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1. INTRODUCTION 


Pulses are unique crops that are a rich source of vegetable protein and have in-built 
mechanism to fix atmospheric nitrogen. India is a major pulse growing country in the world, 
with about 37% of area and 28% of the production. Pulses or food legumes are the most 
versatile and nutritious food available that includes beans, peas, and lentils. They are 
typically low in fats, contain no cholesterol and are high in proteins, potassium, iron and 
magnesium. They also have phytochemicals; a group of compounds that are known to help 
prevent chronic diseases such as cardiovascular diseases and cancers. In addition, they are a 
good source of fiber that can reduce the risk of developing diabetes and help lower blood 
cholesterol levels, thereby reducing the risk of heart diseases. Their by-products provide 
nutritious fodder to livestock. Because of their ability to fix atmospheric nitrogen, they play a 
key role in maintaining soil fertility and ensuring sustainability of production system, 
particularly in low input, small-scale agriculture. ). India is the largest producer, importer and 
consumer of the pulses in the world accounting for 25% of the global production, 1 5 % of the 
trade and 27% consumption, as sizeable population of the country still depends on the 
vegetarian diet to meet its protein requirement. The country produces variety of pulses 
including chickpea (40%), pigeon pea (18%), urd bean (11%), mung bean (9%), lentil 
(8.6%), field pea (5%) and others to the tune of 1 3 to 15 million tones from an area of 22-33 
million hectare with an average yield of 600 to 650 kg/ha. Domestic production of pulses 
after its peak of 14.91 million tones in 2003-04 had declined to 13.13 million tones in 2004- 
OS and to 13.39 million tones in 2005-06 (Ministry of Agriculture) due to adverse climatic 
conditions prevalent in the major production zones. In the current year (2006-07), the 
production is estimated to increase marginally to 14.1 million tones. This still falls short of 
the domestic requirement of 1 7 million tones which is increasing consistently with growing 
population, rising income, value addition and market opportunities. The widening gap in 
demand and supply has led to soaring prices of pulses during past two years. Also, inclusion 
of pulses in future trading and limited availability in the international market has further 
fueled the prices in the prices. 




Molecular Mapping 


A population used for gene mapping is commonly called a mapping population. Mapping 
populations are usually obtained from controlled crosses. Decisions on selection of parents 
and mating design for development of mapping population and the type of markers used 
depend upon the objectives of experiments, availability of markers and the molecular map. 
The parents of mapping populations must have sufficient variation for the traits of interest at 
both the DNA sequence and the phenotype level. The variation at DNA level is essential to 
trace the recombination events. The more DNA sequence variation exists, the easier it is to 
find polymorphic informative makers. When the objective is to search for genes controlling a 
particular trait, genetic variation of trait between parents is important. If the parents are 
greatly different at phenotypic level for a trait, there is a reasonable chance that genetic 
variation exists between the parents, although uncontrolled environmental effects could 
create large phenotypic variation without any genetic basis for the effects. However, lack of 
phenotypic variation between parents does not mean that there is no genetic variation, as 
different sets of genes could result in same phenotype. 

Selection of Parents for Developing Mapping Population 

Selection of parents for developing mapping population is critical to successful map 
construction. Since a map’s economic significance will depend upon marker-trait association, 
as many qualitatively inherited morphological traits as possible should be included in the 
genetic stocks chosen as parents for generating mapping population. Consideration must be 
given to the source of parents (adapted vs 1 exotic) used in developing mapping population. 
Chromosome pairing and recombination rates can be severely disturbed (suppressed) in wide 
crosses and generally yield greatly reduced linkage distances (Albine and Jones, 1987; Zamir 
and Tadmar, 1986). Wide crosses will usually provide segregating populations with a 
relatively large array of polymorphism when compared to progeny segregating in a narrow 
cross (adapted x adapted). To have significant value in crop improvement programme, a map 
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made from a wide cross must be collinear (i.e. order of loci similar) with map constructed 
using adapted parents. 

Types of mapping populations 

Different types of mapping populations that are often sed in linkage mapping are: (i) F2 
population; (ii) F2 derived F3 (F2:F3) populations; (iii) Backcrosses; (iv) Doubled haploids 
(DHs); (v) Recombinant Inbred Lines (RILs); and (vi) Near-isogenic Lines (NILs). The 
characteristic features, merits and demerits of each of these populations are briefly presented 
below: 

F2 population : 

> Produced by selfing or sib mating the individuals in segregating populations 
generated by crossing the selected parents. 

> F2 individuals are products of single meiotic cycle 

> Ratio expected for dominant marker is 3 : 1 and for codominant marker is 1:2:1 

Blackgram or urdbean is one of the important grain legumes of India and Munhgbean yellow 
mosaic virus (MYMV) disease is one of the most devastating diseases. Munhgbean yellow 
mosaic virus (MYMV), a whitefly-transmitted bigeminvirus, causes disease in blackgram. 
This virus is assigned to the genus Begomovirus within the family Geminiviridae. Plant 
resistance is the most common and efficient method of MYMV control. In blackgram, two 
symptom types yellow mosic and necrotic mottle, can be distinguished (Nair and Nene, 
1974). Other crops, including cowpea, develop similar symptoms as a result of infection with 
MYMV. The necrotic mottle is usually associated with resistance (varma et al. 1992). 

Linkage map construction typically requires cosegregation analyses of hundreds of 
Mendel i an loci, most of them molecular markers, using a relatively large number of plants 
from a population in linkage disequilibrium (usually F2, BC1, or similar progenies). The 
mapping effort is considerable, particularly when the objective is to obtain a high-density 
map or to incorporate a large number of functionally meaningful markers, such as those 
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based on expressed sequence tag (EST) sequences, or markers that are particularly suitable 
for breeding applications, such as simple-sequence repeat (SSR) markers, into an existing 
map. A strategy to improve the efficiency of mapping, named selective mapping, was 
proposed by Visionet al. (2000). It consists of a two-step process in which, first, a mapping 
population of usual size (N1460-250) is used to construct a saturated framework map with 
markers placed on it with high precision, and second, new markers are added to this map 
with lower precision using a selected subset of highly informative plants. The final objective 
is to lower the cost of genotyping new markers with a minimal loss of mapping precision. 
The selection of this subset of plants is based on the number and position of recombinational 
crossover sites (or breakpoints) detected with the framework marker data in each plant. The 
breakpoints identified by the ensemble of the selected plants define a set of bins, i.e., 
chromosome fragments bounded by two adjacent breakpoints or by a distal breakpoint and 
the telomere, characteristic of each subset. For a given marker, the joint genotype of the 
selected subset of plants ideally identifies a unique small bin in the genome. The optimal 
subset of a given size would have the maximum possible number of breakpoints evenly 
spaced throughout the genome, resulting in a high number of small bins of uniform size. 
Vision et al. (2000) developed methods and designed a software program (MapPop) to 
facilitate the selection of optimal (or nearly optimal) subsets from mapping populations. 

Diversity analysis 

Among pulse crops, Asiatic grams are assigned to subgenus Ceratotropis (Piper.) Verde, and 
constitute an economically important group of cultivated and wild species for which rich 
diversity found in India. Vigna of the family Fabaceae includes many important pulse- 
yielding taxa. The genus Vigna was originally published by Savi in 1824, which named it 
after Domenico Vigna, Professor of Botany in Pisa (Baudoin and Marechal, 19880ut of 
seven ( V. radiata (mungbean), V. mungo (urdbean), V. angularis (adzuki bean), V. umbellata 
(rice bean), V. aconitifolia (moth bean), V. trilobata and V. trinervia ) domesticated or 
semidomesticated forms, V. radiata and V. mungo are of Indian origin and economically very 
important. These two species are morphologically very similar and only breeding test show 
that these two grams have almost completely separate gene pool. For a long time, V. 
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sublobata was considered as common ancestor of both cultigens, until Arora et al., (1973) 
reported that it includes two distinct wild forms which are not freely cross compatible with 
each other but could be crossed easily with cultigens, one with V. radiata and another with V. 
mungo. Lukoki et al., 1980 reviewed the status of wild Vigna population and suggested that 
two distinct types should be treated as botanical varieties within the respective cultivated 
species namely, V. radiata var. sublobata and V. mungo var. silvestris. However, Babu et al. 
(1985) did not recognize such varietal status for these wild forms and described them as one 
species, viz. V. sublobata chiefly, due to continuous and more or less overlapping 
morphological variations observed in natural populations in Indian subcontinent. Though the 
confusion regarding the taxonomic status of mungbean and urdbean has been settled in an 
undisputed manner, the identity, taxonomy and their relationship with the wild forms remain 
ambiguous till date. Also, due to the threatened habitats of the wild forms and increasing 
realization of wild relatives as potential reservoirs of useful genes, importance of genetic 
diversity analysis of the wild forms cannot be overemphasized. 

Black gram ( Vigna mungo (L.) Hepper is an important crop of the Vigna group. It is also 
known as urdbean, urad dal, urd bean, urd, urid, black matpe bean, black gram, black lentil. It 
has a strategic position in Southeast Asian countries for nutritional security and sustainable 
crop production (Ali and Kumar, 2006). 


Characterization and analysis of genetic diversity are central to effective management of 
plant genetic resources. Characterization provides description of the germplasm while 
analysis of genetic diversity gives information about the range of genetic variation. Proper 
characterization helps in unambiguous discrimination between accessions, detecting 
redundancies and in monitoring genetic change during maintenance. Estimates of genetic 
diversity are useful to select suitable parents for breeding programmes and to develop in-situ 
and on-farm conservation strategies (Karihaloo, 2004). 

The eukaryotic genome contains repeat regions of different lengths and motifs. The repeat 
sequences that have less than six base long central motifs are called microsatellites or SSRs. 
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The term ‘microsatellites’ was coined by Litt and Lutty in 1989. These loci contain tandem 
repeats that vary in the number of repeat units between genotypes and are referred to as 
variable number of tandem repeats (VNTRs) (i.e. a single locus that contains variable 
number of tandem repeats between individuals) or hypervariable regions (HVRs) (i.e. 
numerous loci containing tandem repeats within a genome generating high levels of 
polymorphism between individuals). Microsatellites and minisatellites thus form an ideal 
marker system creating complex banding patterns by simultaneously detecting multiple DNA 
loci. Some of the prominent features of these markers are that they are dominant 
fingerprinting markers and codominant STMS (sequence tagged microsatellites) markers. 
Many alleles exist in a population, the level of heterozygosity is high and they follow 
Mendelian inheritance. New alleles with different numbers of repeating units are generated 
by slippage of the DNA polymerase during DNA replication (Tautz et al., 1986) and unequal 
crossing over (Levinson and Gutman, 1987). Various functional roles have been attributed to 
SSRs: for example, hotspots of recombination, regulation and expression of gene and sex 
determination in various eukaryotes. 

Microsatellite sequences converted into PCR-based markers 

The simple repetitive DNA sequences of SSRs, which are abundantly spread throughout the 
genome of eukaryotes {one SSR per 64.6 kb in monocotyledonous and one per 21.2 kb in 
dicotyledonous species (Wang et al ., 1994)}, provide the basis of PCR based co-dominant 
genetic marker system called Sequence Tagged Microsatellite Sites (STMS). The regions 
flanking the microsatellites are generally conserved among the genotypes of the same 
species. Polymerase Chain Reaction (PCR) primers from the regions flanking the repeat 
regions of genome are used to amplify the SSR containing DNA fragments. The length 
polymorphism is created when PCR products from the different individuals vary in length as 
a result of variation in the number of repeat units in the SSR (Powell et al 1996). This 
method includes DNA polymorphism using specific primers designed from the sequence data 
of a specific locus. Primers complementary to the flanking regions of the simple sequence 
repeat loci yield highly polymorphic amplification products. Polymorphisms appear because 
of variation in the number of tandem repeats (VNTR loci) in a given repeat motif. Tri- and 
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tetra nucleotide microsatellites are more popular for STMS analysis because they present a 
clear banding pattern after PCR and gel electrophoresis. However, dinucleotides are 
generally abundant in genomes and have been used as markers e.g. (CA)n(AG)n and (AT)n. 
The di- and tetranucleotide repeats are present mostly in the non-coding regions of the 
genome, while 57% of trinucleotide repeats are shown to reside in or around the genes. A 
very good relationship between the number of alleles detected and the total number of simple 
repeats within the targeted microsatellite DNA has been observed. Thus larger the repeat 
number in the microsatellite DNA, greater is the number of alleles detected in a large 
population. 


Microsatellite markers are a preferred type of DNA marker used for germplasm analysis and 
varietal identification, marker assisted selection and genome mapping because they are locus 
specific, widely dispersed throughout the genome, highly polymorphic due to variation in the 
repeat units and highly informative because of co-dominant nature (Weber, 1994). 
Microsatellites are reported to be more variable than RFLPs or RAPDs and have been widely 
utilized in plant genomic studies. The major drawback of microsatellites is that they need to 
be isolated de novo from species that are being examined for the first time. 

Specific markers like STMS (sequence-tagged microsatellite markers) ALPs (Amplicon 
length polymorphisms) or STS markers have proved to be extremely valuable in the analysis 
of gene pool variation of crops during the process of cultivar development, and classification 
of germplasm (Yang et aL, 1994). These markers are extremely sensitive and can detect 
allelic variability during cultivar development (Yang et al., 1994). STS markers specific to 
chloroplast or mitochondrial DNA have been useful in providing seed and pollen specific 
markers which can be utilized for the detection of length variation at multiple physically 
linked sites and may be used to provide haplotype data and thus genotypically unique 
individual plants (Wu et al., 1 998). Also a comparison of patterns of variability detected with 
biparentally (nuclear) and uniparentally (organellar) transmitted markers can provide 
complementary information to population and evolutionary biologists. Excellent examples of 
this are the Poly A mononucleotide repeats in maize (Powell et al., 1995), Poly (TA/AT) 
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dinucleotide repeats found in liverworts, maize, pea and nonphotosynthetic green plant 
Epifagus virginiana, and a total of 500 chloroplast SSRs identified with repeat motifs greater 
than 10 repeat units in rice, tobacco, black pine, liverwort and maize (Powell et al., 1996). 
Though all these marker types provide valuable information regarding the evolution and 
phylogeny of various species being studied in any given set of samples, the trend is now 
shifting towards the use of ESTs (expressed sequence tags) for such analysis. This may be so, 
because in such studies, one actually looks at the evolution of functional genes (Mason- 
Gamer et al. 1998 and Deshpande et al. 1998). With these views in mind the present treatise 
included the following precise objectives. 

Objectives: 

i) Elucidation of genetics of resistance to YMV diseases in urdbean 

ii) Molecular mapping of genes controlling resistance to YMV. 

iii) Studying genetic relationships among important Vigna cultivars and wild 
species using molecular markers. 

iv) Study of genetic variability in different gene pools based on molecular gene 
sequence. 
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REVIEW OF LITERATURE 


The genus Vigna is an important legume taxon, widely distributed in tropical and 
subtropical regions of both hemispheres. It contains about 1 50 species, divided into seven 
subgenera. Of these, the subgenus Ceratotropis, mainly comprising species found in Asia, 
has been differentiated into morphologically homogenous group with specialized and 
complex floral morphology (Marechal et al, 1980). Vigna mungo (L.) Hepper (urdbean) and 
V. radiata (L) Wilczek (mung bean) are the two most important pulse yielding taxa 
belonging to subgenus Ceratotropis of genus Vigna. The Indian center is extremely rich and 
possesses considerable landrace diversity of these species and several allied ancestral wild 
taxa. V. radiata and V. mungo are morphologically similar, however minor differences exist 
in morphology and important plant characters. Variations have also been reported in 
karyotype, DN A content, proteins, amino acid composition and species hybridity. 

2.1 Taxonomy And Species Relationship 

The taxonomic and nomenclatural confusion occurred due to incorrect nomenclature by 
Linneaeus, who named urd bean as Phaseolus mungo and mung bean as P. radiatus. 
Verdcourt, 1970 reviewed the genus and on the basis of morphological and biochemical 
evidences transferred Asian species of Phaseolus to genus Vigna. Later, Phaseolus mungo 
was named as V. mungo and P. radiatus as V. radiata. Also, it was considered that both these 
species were evolved from common ancestor P. sublobatus, which afterwards named as V. 
sublobata. Also, further misunderstanding occurred due to wrong cytogenitical relationships 
based on cytology study (De & Krishnan, 1966) and scanning electron microscopy of seed 
coat pattern (Jain & Mehra, 1980). The latter author contented that V. sublobata was not 
likely the progenitor of mungbean (V. radiata) but instead the progenitor of urdbean. 
Morphological and biochemical studies by Arora et al, 1973 on wild populations of Vigna 
provide the first constructive evidence about evolutionary relationships among V. radiata- 
mungo-sublobata-silvestris complex. They identified two categories, which respectively gave 
rise to V. radiata and V. mungo. Lukoki et al, 1 980 accepted specific distinction between two 
forms and their relation as wild ancestor to the cultivated species. They noted that Plx 416 
was similar to V. mungo biochemically in possessing glutamyl-methionine and its sulphoxide 
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and morphologically possess narrower stipules, brighter yellow flowers, 6-8 ovules per 
ovary, erect pods with long white hairs (or glabrescent) and seeds with a raised arillate hilum. 
Therefore, they recognized these plants as the wild ancestor of black gram. It was described 
as a new taxon, namely V. m. var. silvestris Lukoki, Marechal & Otoul. Lukoki et al (1980), 
also described Plx 274 similar to V. radiata var. radiata in possessing a-glutamyl-s- 
methylcysteine and its sulphoxide in its seeds and have broader stipules, paler yellow 
flowers, 10-14 ovules per ovary, spreading pods with short brown hairs and seeds with a flat 
non-arrillate hilum. Chandel (1984) and Miyazaki (1984) also supported the view that 
mungbean and urdbean were domesticated from V. radiata var. sublobata and V. m. var. 
silvestris respectively on the basis of their biochemical and morphological studies. 

2.2 Origin 

V. radiata, V. mungo and their proposed wild forms are originated in India. Archaeological 
records for these species were not found anywhere outside India (Kajale, 1974). Charred 
grains of mungbean were recorded from Navadatoli (1500-1000 BC), Satavahana Bhokardan 
(200 BC- 200 AD) and Neolithic Chirand (1800-1300 BC) in Bihar. Carbonized grains of V. 
mungo were discovered from Chalcolithic Navadatoli (1500-1000 BC), Neolithic Piyampalli 
(1400 BC) in Tamilnadu and Indo Roman Nevasa (50 BC-200 AD) in Maharashtra. Also, 
charred grains of wild and cultivated species occurring in stratified layers (Kajale, 1977) 
provide evidences of various stages of evolution from wild to domesticated forms. 

Kajale, 1 974 pointed out that phytogeographic distribution of ancient grains might not 
correspond with the present day distribution of wild ancestral forms, or even cultigens. 
However, Chandel (1984) showed that archaeological evidences correspond very well with 
the distribution and co-occurrence of present day putative progenitors of mungbean and 
urdbean. The wild ancestral form of urdbean, V. mungo var. silvestris corresponded with the 
archaeological remains at Navadatoli and Bhokardan, while charred grains of mungbean 
found at Chalcolithic Navadotoli and Neolithic Chirand could be correlated with the 
occurrence of V. radiata var. sublobata . Also, there is mention of these species in Vedic text 
such as Kautilya ‘Arthshastra’ and in ‘Charak Samhita’. These evidences provide strong 
support for the origin and domestication of mung bean and urd bean in Indian gene center. 
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2.3 Ecology, Distribution And Domestication 


V. radiata and V. mungo are cultivated throughout Asia but are characteristically adapted to 
particular and different agro climatic conditions. V. radiata is widely distributed, relatively 
drought tolerant and well adapted to a range of soil conditions including light soils and 
limited irrigation. It is cultivated throughout the southern half of Asia, especially India, 
Bangladesh, Pakistan, Sri Lanka, Thailand, Laos, Cambodia, Vietnam and relatively dry 
parts of Java, Eastern Malaysia, South China and Central Asia. V. mungo on the other hand 
require wetter conditions and thrive in relatively heavier soils. It is however unsuitable for 
the wet tropics and grows best, where rainfall is less than 1000 mm per annum. As a crop it is 
limited to parts of Pakistan, India, Sri Lanka and Burma, but it has recently been taken up in 
parts of Southeast Asia and elsewhere in the tropics, by Indian immigrants. V. radiata var. 
sublobata is the most widely distributed wild species of the subgenus Ceratotropis. It grows 
along the forest margins and grasslands and in disturbed habitats such as roadsides. It is 
found across equatorial Africa as far west as Benin, also Tanzania, Madagascar; Middle East- 
Oman; South Asia-India, Bangladesh, Sri Lanka; South East Asia - Myanmar, Thailand, 
Vietnam, the Philippines, Indonesia; East Asia-China including Taiwan; Papua New Guinea 
and Australia. V. mungo var. silvestris grows in disturbed habitats such as, roadsides. It is 
found only in India, Myanmar and Thailand. 

Domestication is a multistage process, for pulses it includes increase in seed size, 
reduced pod shattering and importantly loss of germination inhibition (Plitman and Kislev, 
1989; Smartt, 1990; Zohary and Hopf, 2000). Archaeology together with the modem day 
distribution of wild progenitors suggests three distinct centers domestication in India, 
including South India, Gujarat and the Ganges plain (Fuller 2002, 2006). These 
domestication events were associated with mid-Holocene, between 4000 and 2000 BC, and 
occurred amongst local populations of hunter-gatherers in monsoonal India, although some 
influences from the west due to spread of livestock cannot be ruled out. The V. radiata was 
probably brought under cultivation in the Northwest, probably in the Himalayan foothills of 
Punjab region, and in the far south, as early finds in both regions are widely separated. Early 
finds of urdbean come Gujarat and the Northern Peninsula where wild population of these 


11 



species persists (V. mungo var. silvestris) without association with wild mungbean (V. 
radiata var. sublobata) (Fuller and Harvey, 2006; Tomooka et al, 2003). Early finds in South 
India of V. radiata occur in the driest Savannah environments of Peninsula together with 
large quantities of Horsegram ( Macrotyloma uniflorum ) and small millets ( Brachiaria 
ramosa, Setaria verticillata ), and in some cases associated with introduced domesticates 
(Fuller et al, 2004). This association with other crops, there large quantity and the difference 
from the habitat of the wild progenitor all imply cultivation of these early mungbeans (Fuller 
and Korisetter, 2004). However, wild form of mungbean is widely distributed in Africa, Asia 
and Australia, so domestication more than once cannot be ruled out. 

Tomooka et al, (1992) based on seed protein variation in mungbean landraces from Asia, 
identified regions of protein diversity and proposed two paths of dissemination. Protein type 
diversity is greatest in West Asia (Afghani stan-Iran-Iraq) rather than India. From this region 
mungbean was proposed to have moved east by two routes one following a Southern route 
through India to Southeast Asia. The other pathway followed a northern route along the Silk 
Road or possibly through India to China. 

2.4 Cytology, Crossability And Species Hybridization 

Several cytological and cytogenetical studies have been carried out to determine 
chromosome number (Bose, 1939; Sen and Ghosh, 1961; Coyne, 1964) and chromosome 
homologies (De and Krishnan, 1966, 1968a, 1968b) in genus Vigna. Subsequently series of 
investigations were made on Cytology, Crossability and Species hybridization (De and 
Krishnan, 1968b; Dana, 1966a, b, c, d; Biswas and Dana, 1976a, b; Machando et al, 1982). 
Both mung and urd bean are diploid with chromosome number 2n=2x=22. The investigation 
related to the organization and evolution of genome of Vigna species appeared difficult due 
to size, large chromosome number and uniformity in shape and size within and between the 
complements. 

Panda et al 1991, studied quantitative DNA variations between and within chromosome 
complements of Vigna species. The different accessions of V. radiata did not exhibit 
significant differences in 2C DNA content. The DNA amount in V. radiata narrowly ranged 
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from 2.63 to 2.70 pg (mean 2.672 pg), whereas, V. mungo possessed 2.83 pg. A series of 
investigations were conducted involving both cultivated and wild progenitors (Rao and 
Chandel, 1991; Chandel et al, 1991; Sehgal and Chandel, 1991; Audilakshmi and Chandel, 
1990). Chromosome number, morphology, nucleolar organization, karyotype 
symmetry/asymmetry of V. radiata, V. radiata var. sublobata, V. mungo and V. mungo var. 
silvestris were studied. The study revealed significant differences in the karyotypes, both in 
between and within the cultivated as well as wild progenitor forms. Variations occurred in 
the number of V, L and J chromosomes. Rao and Chandel, 1991 confirmed in their studies 
that species/forms of Vigna occurring in Indian gene center have evolved directly from their 
ancestral forms (Chandel, 1984; Chandel et al, 1984) without undergoing any numerical 
changes in their chromosome complements. The investigation supported the view that basic 
chromosome number of mungbean and urdbean is x=T 1. Beside there is overall stability in 
chromosome morphology and symmetry in cultigens and wild alien species, morphological 
diversity could not be directly correlated with the karyological differentiation. It was also 
concluded that V. radiata var. sublobata is surely the wild progenitor of mung bean and 
definitely not of urd bean (V. mungo ) as reported earlier by Sharma et al, 1980 and Babu et 
al, 1988. 

Crosses performed between mung and urd bean (Dana, 1966a; Sen and Bhowal, 
1960) showed presence of sterile hybrids. In these crosses translocations, inversions, 
deficiencies and duplications in chromosomes were suggested as the possible reasons for 
hybrid inviability. Similarly, when V. radiata var. sublobata was crossed with V. radiata it 
produced fertile hybrids, while when crossed with V. mungo it resulted in only few shriveled 
seeds (Karmakar and Dana, 1987). They attributed it to cytogenetic differentiation between 
V. radiata and V. mungo that appears to have occurred during process of natural evolution 
resulting in the speciation of these two very similar taxa. Also, Karmakar and Dana, 1987 
reported that wild collection of V. radiata var. sublobata (Morphotype 2-Jabalpur) when 
crossed with V. radiata as female parent, they found that V. radiata var. sublobata was not 
separated from mung bean by any isolating barriers, while it is separated from urd bean by 
multiple types of isolating mechanisms. 
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Comprehensive interspecific hybridization studies were undertaken between V. radiata and 
its allied wild prototype species V. radiata var. sublobata, V. mungo and its wild ancestral 
form V. mungo var. silvestris (Audilakshmi and Chandel, 1990). Different species/accessions 
were involved in straight as well as in different combinations including reciprocal crosses. 
The close relationship between wild progenitor species and their cultigens was clearly 
established with the production of successful hybrids and normal segregations in subsequent 
generations. However, V. radiata var. sublobata exhibited only partial compatibility with V. 
mungo , suggesting distinctness. These results were confirmed through meiotic studies as 
well. 

2.5 Biochemical Analysis 

The electrophoresis of seed proteins showed marked differences in mungbean and urdbean 
employing gel electrophoresis technique (Mihardja et al, 1974). Lukoki et al 1980, showed 
presence of glutamyl 1-s-methylcysteine in mungbean and its wild progenitor V. radiata var. 
sublobata while its absence was shown in V. mungo and its putative wild progenitor V. 
mungo var. silvestris. They also showed absence of glutamyl methionine and its sulfoxide in 
V. radiata and V. radiata var. sublobata complex, while above-mentioned amino acids were 
present in V. mungo and V. mungo var. silvestris. Chandel (1980) carried out comprehensive 
biosystematics studies involving nine cultivated and wild species of Vigna. The author 
carried out comparative morphology, scanning electron microscopy of seed coat cellular 
structure and hilum structure, chromatography of leaf phenolics and gel electrophoresis of 
seed proteins. The chromatographic studies revealed a total of 62 phenolic spots in all nine 
species. V. radiata exhibited presence of 17 spots while V. mungo possessed 21 spots. V. 
radiata var. sublobata exhibited similar phenolic compounds and their position, and Rf value 
corresponds with mungbean. Also, wild ancestral form V. mungo var. silvestris showed 
remarkable similarity with its putative domesticated form V. mungo. 

Babu et al, 1988 assessed the seed protein and amino acid composition of V. radiata 
var. sublobata and two cultivated forms V. mungo and V. radiata. They reported high range 
of variation for all the amino acids indicating broad genetic base. They suggested the 
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usefulness of wild populations of V. radiata var. sublobata in the nutritional upgrading of 
mung and urd bean. The amino acid profiles were population specific and wild population of 
V. radiata var. sublobata were shown to differ contrastly in cysteine content with V. radiata 
as well as V. mungo. V. mungo contained high methionine content, which could be due to Y- 
glutamylmethonine and its sulphoxide while V. radiata possessed cysteine and its 
sulphoxide. These results were similar to that reported by Lukoki et al, 1980 and Otoul et al, 
1975. However, Babu et al continued to stress upon that V. mungo and V. radiata var. 
sublobata are closely related. It may be due to the fact that the population of V. sublobata 
they have analyzed could be that of V. mungo var. silvestris as it was collected from Western 
Ghats and therefore resembled more to V. mungo. 

Some indirect evidences from nodulation studies also indicate the origin of V. mungo 
and V. radiata from V. mungo var. silvestris and V. radiata var. sublobata respectively 
(Kavimandan and Chandel, 1988). Several strains of root nodule bacteria were isolated form 
V. mungo, V. radiata, V. radiata var. sublobata and V. mungo var. silvestris. Wild progenitor 
of V. radiata and V. mungo; V. radiata var. sublobata and V. mungo var. silvestris 
respectively exhibited distinct differences for bacterial specificity resulting in differential 
nodulation and the nitrogen accumulation. 

2.6 Genomic Characters, Photosynthetic Activity, Seed Protein Content And Amino 
Acid Levels 

Genomic characters such as the variations in the 2C DNA content are related to duration of 
cell cycle, minimum generation time, latitudinal distribution of agricultural crops and nuclear 
instability of hybrids (Benett, 1976, 1977, 1982; Smith and Benett, 1975; Rees et al, 1982). 
The 2C DNA amount was low in wild forms in comparisons to cultivars (Ignacimuthu and 
Babu, 1988b; Lakhanpaul, 1989). This suggests that evolution of 2C DNA in genus Vigna is 
towards an increase in 2C DNA content. This increase is probably associated with 
domestication and selection. The wild forms have to adapt to fluctuations in environmental 
conditions; therefore the lower 2C DNA content enable the plant populations to have shorter 
mitotic and meiotic cycles. On the other hand, cultivars have been evolved as a result of 
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directional selection under stable environment; as such the longer duration of mitotic and 
meiotic cycles may not be a major constraint in the reproductive efficiency of the agricultural 
populations (Nagato et al, 1981; Resslar et al, 1981). Also, induced chromosomal aberrations 
in wild relatives are low in comparison to cultivars (Kallo, 1972; Wolff, 1981; Gupta and 
Roy, 1985; Ignacimuthu and Babu, 1989). This implies that the genomes of the cultivars are 
rapidly differentiated because of larger number of macro mutations leading to higher rates of 
sequence turnover (Flavell, 1982; Rees et al, 1982; Rees, 1984). Domestication and selection 
might have further aided in the rapid differentiation of genomes between cultigens and 
divergence of the cultigens from wild relatives. 

Experiments on photosynthetic activity by Ignacimuthu and Babu (1987, 1989), 
showed that wild relatives have high photosynthetic rate and RuBPCase activity at vegetative 
stage as compared to the cultivated forms. This might be due to the high rate of 
photosynthesis per unit leaf area in wild forms in comparison to the cultigens (Evans, 1983). 
The cultivars have shown higher photosynthetic activity the reproductive stages due to their 
larger sink size than the wild relatives (Lorimer, 1981; Me Cree, 1983). The higher 
photosynthetic rate at the vegetative stage in wild forms was due to smaller leaf area and 
higher dry weight of leaf (Ignacimuthu, 1985). The increase in leaf area coupled with the 
increase in size and number of harvestable organs during the process of selection and 
domestication led to increase the photosynthetic rate of cultigens, but the Carbon dioxide 
Exchange Rate (CER) per unit leaf area is low because of larger leaf size (Gifford and Evans, 
1981; Yamauchi and Yoshida, 1985). The wild species have productivity promoting gene 
assemblies and hence of immense value in the development of high yielding varieties 
(Ignacimuthu and Babu, 1987). 

Studies on seed proteins revealed that wild forms are as good as cultivars, 
nutritionally (Ignacimuthu and Babu, 1984, 1987; Babu et al, 1988). Also the amino acids 
levels of the wild forms are higher than the FAO standards and are even higher than the 
cultivated forms (Ignacimuthu and Babu, 1987, 1989; Babu et al, 1988). This genetic 
variability found in essential amino acids of seed proteins among different population could 
be used in cross breeding programmes for the nutritional upgrading of modem cultivars 
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(Gepts and Bliss, 1984). Therefore, the wild populations are nutritionally as good as or even 
better than the cultigens grown for the human consumption or livestock feed. 

2.7 Wild Relatives As Genetic Resources 

Fuzi and Miyazaki (1987) reported an accession (TC1966) of V. radiata var. sublobata that 
showed perfect resistance against azukibean weevil ( Callasabruchus chinensis ). The 
resistance was found to control by a single dominant (Kitamura et al, 1988). Fuzi et al (1989) 
further found that TCI 966 is completely resistant against C. maculates, C. phaseoli and 
Zabrotes subfasciatus. Tomooka et al (1992) developed a bruchid resistant mungbean line in 
Thailand by using TCI 966 as a gene source. In addition to bruchid resistance, high resistance 
to yellow mosaic virus (Singh and Ahuja, 1977), high methionine content in seeds (Babu et 
al, 1988), higher photosynthetic activity and tolerance to drought ((Ignacimuthu and Babu, 
1987), higher tolerance to saline and alkaline soils (Lawn et al, 1988) have been reported for 
V. radiata var. sublobata. In comparison to V. radiata var. sublobata there are very few 
studies of evaluation of V. mungo var. silvestris as genetic resource. However, this variety is 
cross compatible with V. mungo (Miyazaki et al, 1984). Therefore, further collection and 
evaluation V. mungo var. silvestris is necessary for genetic upgradation of V. mungo. 

The results of above mentioned studies conducted by several independent workers 
and support from taxonomic, cytological/cytogenetical and species hybridization studies as 
well as evidences from biochemical assays, SEM studies and archeological sources indicate 
that V. radiata and V. mungo are two distinct species. V. radiata and V. mungo have evolved 
directly from their proposed wild forms V. radiata var. sublobata and V. mungo var. silvestris 
respectively in Indian gene center. 


2.8 DNA Based Molecular Markers In Plants 

With the advent of DNA based molecular markers over last two decades, the entire scenario 
of biological sciences has been revolutionalized. Molecular markers have acted as versatile 
tools and have carved their own niche in various fields like taxonomy, physiology, genetic 
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engineering etc. They are no longer merely looked upon as DNA fingerprinting markers in 

variability studies or as forensic tools. Ever since their development, they are constantly 

being modified, to enhance their utility and to bring about automation in the process of 

< 

genome analysis. The discovery of polymerase chain reaction (PCR) was landmark in this 
effort and proved to be a unique process that opened up a multitude of new possibilities in 
the field of molecular biology. This has facilitated the development of marker based gene 
tags, map-based cloning of agronomically important genes (Peng et al, 1999), variability 
studies (Fahima et al, 2002; Thiang et al, 2002) and phylogenetic analysis (Matsuoka et al, 
2002a); thus giving new dimensions to concerted efforts of breeding and marker-aided 
selection for development of newer and better varieties in short time span. These DNA 
markers offer several advantages over traditional phenotypic markers, as they provide data 
that can be analyzed objectively. Thus it can be discerned that DNA based markers has 
brought in the much needed amalgamation of classical breeding and modem biotechnological 
approaches. 

Genetic polymorphism is classically defined as the simultaneous occurrence of a trait 
in the same population of two or more discontinuous variants or genotypes. Although DNA 
sequencing is a straightforward approach for identifying variations at a locus, but until entire 
genomes are sequenced, geneticists will continue to rely on genetic markers as the means of 
characterizing the genotypic variation of individuals. A wide variety of techniques have, 
therefore, been developed in the past few years for visualizing DNA sequence polymorphism 
and are generally classified as hybridization-based markers and PCR-based markers. In the 
former, DNA profiles are visualized by hybridizing the restriction enzyme digested DNA to a 
labeled probe, which is a DNA fragment of known origin or sequence (e.g. RFLP). PCR- 
based markers involve in vitro amplification of particular DNA sequence or loci, with the 
help of specifically or arbitrarily chosen oligonucleotide sequences 9primers) and a 
thermostable DNA polymerase enzyme. The amplified fragments are separated 
electrophoretically and banding patterns are detected by different methods such as staining 
and autoradiography. 
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2.8.1 Random Amplified Polymorphic DNA 

Random Amplified Polymorphic DNA (RAPD) involves the use of single arbitrary primer in 
a PCR reaction, resulting in amplification of many discrete DNA products revealing 
nucleotide sequence polymorphisms. The technique was developed independently by two 
different laboratories (Williams et. al., 1990; Welsh and McClelland, 1990) and called as 
RAPD and AP-PCR (Arbitrary primed PCR) respectively. In this reaction, a single species of 
primer binds to the genomic DNA at two different sites on opposite strands of the DNA 
template. If these priming sites are within an amplifiable distance of each other, a discrete 
DNA product is produced through thermocylic amplification. The polymorphisms between 
individuals result from sequence differences in one or both of the primer binding sites, and 
are visible as the presence or absence of a particular RAPD band. Such polymorphisms thus 
behave as dominant genetic markers. 

The RAPD reaction mixture comprises of a target DNA sequence that is 
exponentially amplified with the help of arbitrary primers, a thermostable DNA polymerase, 
deoxy nucleotide triphosphates, magnesium and reaction buffer. The reaction involves 
repeated cycles, each consisting of a denaturation, a primer annealing and an elongation step. 
In the first step the DNA is made single stranded by raising the temperature to 94°C 
(denaturation). In the second step, lowering of the temperature to about 32° to 45° C results in 
annealing of the primer to their target sequences on the template DNA (annealing step). In 
the third cycle, temperature is chosen where the activity of the thermostable Taq DNA 
polymerase is optimal, i.e., usually 72(C. The polymerase now extends the 3’ ends of the 
DNA - primer hybrids towards the other primer - binding site. Since this happens at both 
primer-annealing sites on both the DNA strands, the target fragment is completely replicated. 
Repeating these three step cycles 40 to 50 times results in the exponential amplification of 
the target between the 5' ends of the two primer binding sites. Amplification products are 
separated by gel electrophoresis and visualized by ethidium bromide staining. 

The RAPD technology has provided a quick and efficient screen for DNA sequence 
based polymorphisms at a very large number of loci. The major advantage is that no DNA 
sequence information of target genome is required. Sets of short primers (usually 10 mers) 
suitable for RAPD amplification are available commercially or can easily be synthesized and. 


19 



apart from a thermocycler and an agarose gel assembly, no special equipment is required. 
When compared with RFLP, RAPD requires small amount of DNA and often-crude miniprep 
procedures yield DNA of sufficient quantity and quality. It does not require radioactivity or 
special staining procedures to visualize polymorphisms. Automation is possible at all stages 
from DNA extraction to data collection and analysis. However, one of the inevitable trade- 
offs with the technique is that amplification is performed under conditions of low stringency. 
Consequently, some of the products result from weak complexes between primer and 
template, and this can result in poor reproducibility for some primers and bands. Another 
problem reported for RAPD analyses is a low incidence of non-inherited bands. While great 
majority of RAPD bands are known to be inherited as Mendelian markers, care is needed 
when drawing conclusions based on a small number of band differences. A limitation of 
RAPD that cannot be overcome is band dominance. This means that heterozygotes are rarely 
detected and there are usually only two states for a polymorphism: present and absent. 
Consequently, RAPD and other multilocus profiles provide less genetic data than profiles for 
single locus codominant markers such as STS. 

RAPD markers have been used for the development of genetic maps, for tagging 
useful traits and in population genetics. One of the first practical uses of RAPD markers was 
in the creation of high-density genetic maps. For mapping with RAPD markers it is necessary 
to use backcross or recombinant population, haploid or gametophytic tissue, or alternatively 
in a F2 population where only RAPD markers amplified from a single parent are mapped. 
Genetic simulations show that dominant markers linked in coupling are as efficient for 
mapping as codominant markers on a per gamete basis. 

RAPD has been used to identify molecular markers that lie within regions of a 
genome introgressed during the development of near isogenic lines. In tomato, RAPD 
markers specific to chromosome 6 were identified by screening Lycopersicon esculentum 
substitution line (Klein - Lankhorst et. al. 1991). RAPD markers linked to the Pto locus in 
tomato were identified after screening two near isogenic lines (Martin et. al., 1991). 

RAPD markers can also be used to screen for markers linked to specific regions of 
the genome by the bulk segregant analysis method. This method uses two-bulked DNA 
samples gathered from individuals segregating in a single population. Each bulk is composed 
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of several individuals that differ for a specific phenotype or genotype, or individuals at either 
extreme of a segregating population. Random primers amplify loci from each pool and 
identify polymoiphisms linked to the trait of interest. 

The area of research that has shown the most growth with respect to the use of RAPD 
technology is that of population genetics. RAPD markers have been used to create 
fingerprints for the study of individual identity and taxonomic relationship in both eukaryotic 
and prokaryotic organisms. RAPD markers are being used effectively to assess the amount of 
genetic diversity in germplasm collections. Using only 25 different decamer oligonucleotide 
primers, Kresovich et al. (1992) collected information on 140 different polymorphic 
characters in a "test array" of individuals representing B. oleracea and B. rapa. They showed 
the utility of the assay for discriminating between individuals in a germplasm collection and 
the ability to distinguish between closely related individuals was simply a function of the 
number of RAPD bands that were observed. 

2.8.2 Inter Simple Sequence Repeats 

Inter Simple Sequence Repeats (ISSR) technique is a PCR-based multi-locus marker system 
that employs oligonucleotide primers homologous to SSR sequences such as (GATA)n to 
amplify mainly the inter-SSR regions. Amplification products are only obtained only if SSRs 
in opposite orientation are found in a PCR-able distance, with flanking sequences matching 
the oligonucleotides. ISSR markers are multilocus unmapped dominant markers. The various 
synonyms used for ISSR techniques are SSR-anchored PCR; anchored SSR-PCR (ASSR- 
PCR); microsatellite primed PCR (MP-PCR); anchored microsatellite primed PCR (AMP- 
PCR); single primer amplification reaction (SPAR); randomly amplified microsatellites 
(RAM) and randomly amplified microsatellites polymorphism (RAMP). The microsatellite 
repeats used, as primers can be dinucleotide, trinucleotide, tetranucleotide and 
pentanueleotide. The primer used can either unanchored (Gupta et al, 1994; Meyer et al, 
1993; Wu et al, 1994) or more usually anchored at 3' or 5' end with 1 to 4 degenerate bases 
extend into flanking sequences (Zuetkiewicz et al, 1994). When unanchored i.e. only SSRs 
are used as primers, the primer tends to slip within the repeat units during amplification 
leading to smears instead of clear bands. Extending the primer (anchoring) with 1 to 4 
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degenerate nucleotides the 3' end or 5' end assures annealing only to the ends of the micro 
satellite template DNA thus obviating internal priming and smear formation. Secondly, the 
anchor allows only a subset of microsatellite to serve as priming sites. When 5' anchored 
primers are used, the amplified products include the microsatellite sequences and their length 
variations across a genome and therefore give more number of bands and higher degree of 
polymorphism. Usually dinucleotide repeats, anchored either at 3' or 5 '-end reveal high 
polymorphism (Blair et al., 1999; Josh! et al., 2000; Nagaoka and Oghihara, 1997). The 
primer anchored at 3 '-end give clear banding pattern as compared to those anchored at 5' 
ends (Tsumura et al., 1996; Blair et al, 1999; Nagaoka and Oghihara, 1997). Since the primer 
is a SSRs motif the frequency and distribution of the microsatellite repeat motifs in different 
species also influence generation of bands. 

ISSR technique is nearly identical to RAPD technique except that ISSR primer 
sequences are designed from microsatellite regions. High stringency of amplification in the 
form of long primers (16-25 bp) and elevated annealing temperature (45-60° C) makes ISSR 
more reproducible than RAPD markers. ISSRs segregate mostly as dominant markers 
following simple Mendelian inheritance (Tsumura et al., 1996; Gupta et al., 1994; Wang et 
al., 1998; Ratnaparkhe et al, 1998). However they have also been shown to segregate as co- 
dominant markers in some cases thus enabling distinction between homozygotes and 
heterozygotes (Wu et al., 1994; Akagi et al., 1996; Wang et al, 1998; Shankar and Moore, 
2001 ). 

ISSR finds application in genomic fingerprinting (Charters and Wilkinson, 2000; 
Prevost and Wilkinson, 1999; Wolff et al, 1995), genetic diversity and phylogenetic analysis 
(Virk et al, 2000; Qian et al, 2001; Yang et al, 1996; Wang et al, 1998), genome mapping 
(Backer et al, 1995; Kojima et al, 1998; Arcade et al, 2000), gene tagging and marker 
assisted selection (Akagi et al, 1996; Levin et al, 2000; Hantula et al., 1996), determining 
SSR motif frequency (Blair et al, 1999; Nagaoka and Oghihara, 1997) and studies on natural 
population/speeiation (Qian et al, 2001; Wolff et al, 1995). 
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2.8.3 Amplified fragment length polymorphisms 

Amplified fragment length polymorphisms (AFLPs) (Vos et al., 1995) are Polymerase Chain 
Reaction (PCR) based markers for DNA fingerprinting and rapid screening of diversity. It 
combines reliability of RFLP with the power of PCR. The technique is based on the selective 
PCR amplification of subsets of restriction fragments from a total digest of genomic DNA. It 
comprises of the restriction of DNA with two restriction enzymes preferably a hexa cutter 
and a tetra cutter, the ligation of double stranded (ds) oligonucleotide adapter to the ends of 
the restriction fragments, the amplification of a subset of the restriction fragment using two 
primers complementary to the adapter and a restriction site sequences, and extended at their 
3’ ends by ‘selective’ nucleotides, gel electrophoresis of the amplified restriction fragments 
on denaturing polyacrylamide gels (sequence gels) or on capillary based gels and the 
visualization of the DNA fingerprints by means of autoradiography, phosphoimaging or other 
methods. 

The Amplified fragment length pattern generated by means of the AFLP technology 
are known as AFLP fingerprint. The AFLP fingerprints are rich source of restriction 
fragment polymorphisms, termed AFLP markers. The frequency with which AFLP markers 
are detected depends on the level of sequence polymorphism between the tested DNA 
samples. 

Single nucleotide changes will be detected by AFLP when an actual restriction site is 
affected, or nucleotides adjacent to the restriction sites are affected, which cause the AFLP 
primers to mispair at the 3’ end thus preventing amplification. The selective nucleotides are 
not an exact match of the sequence next to the restriction sites. In addition, deletions, 
insertions and rearrangements affecting the presence or size of restriction fragments will 
result in polymorphism detected by AFLPs. 

The AFLP is the random amplification technique, which in contrast to most other 
random amplification techniques makes use of stringent PCR conditions. The primers are 
generally 17-21 nucleotide in length and anneal perfectly to their target sequences i.e. the 
adapter, restriction site and small number of nucleotides adjacent to the restriction sites. This 
renders AFLP a very reliable and robust technique, which is unaffected by small variations in 


23 



amplification parameters (e.g. thermal cyclers, template concentration, PCR cycle profile). 
The high marker densities that can be obtained with AFLP are an essential characteristic of 
the technology: a typical AFLP fingerprint contains between 50 to 100 amplified fragments, 
of which upto 80% may serve as genetic markers. Moreover, it requires no sequence 
information or probe collections prior to the generation of AFLP fingerprints. AFLP markers 
usually exhibit Mendelian inheritance, indicating that they are unique DNA fingerprints. 

The AFLP markers are used in genetic studies, such as: Biodiversity studies, in 
analysis of germplasm collection, in the genotyping of individuals and genetic distance 
analysis and in the identification of closely linked DNA markers. AFLP markers are also 
used in the construction of genetic DNA marker maps and in the construction of physical 
maps using genomic clones such as YACs and BACs. The use of AFLP in precision mapping 
of genes and the subsequent isolation of genes is also well documented. 

The AFLP is relatively fast (samples can be processed on automated thermal cycler 
and DNA sequencer) and assays the entire genome for polymorphic marker. A relatively 
small amount of genomic DNA (0.05 to 0.5 pg) depending upon the size of genome is 
required. AFLPs provide 10 to 100 times more markers and are more sensitive than other 
fingerprinting techniques. AFLPs are highly reproducible, no prior information of the 
genome is required and also no other taxon specific primer sets are required. Primers are 
commercially available that work for most of the crops. AFLPs are reliable and informative 
multilocus probes (Winfeild et ah, 1998), and provides high levels of resolution to allow 
delineation of complex genetic structures (Powell et al., 1996). A high multiplex ratio (no. of 
different genetic loci analyzed simultaneously per experiment) is possible since AFLP 
markers are distributed across the genome. 

2.8.4 Microsatellites 

Microsatellite are tandem repeats of short (2-5 bp) DNA sequences, the most abundant being 
dinucleotide repeats and represent a significant portion of higher eukaryotic genome (Litt & 
Luty, 1989). These types of reiterated sequences have been termed as differently as Simple 
Sequence Repeats (SSRs, Jacob et al, 1991), Microsatellite (Litt & Luty, 1989) or Short 
Tandem Repeats (Edward et al, 1991). The number of times the unit is repeated in a given 
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microsatellite can be highly variable, a characteristic that makes them useful as genetic 
markers. Majority of microsatellite occur in gene introns or other non-coding regions of the 
genome; thus variation in the number of repeats has no consequence on gene function. 
Generally microsatellite itself does not cause any phenotype, but rather is used as marker to 
identify a specific chromosome or locus. When being used as marker, the specific number of 
repeats in a given microsatellite is not important, but rather the difference in the number of 
repeats between the alleles. The variation in number of repeats affects the overall length of 
microsatellite, a characteristic readily measured by laboratory techniques. This class of 
genetic polymorphism is commonly used for diversity analysis, mapping, linkage analysis 
and tracing inheritance pattern. 

The DNA sequences flanking microsatellite are generally conserved allowing the 
selection of PCR primers that will amplify intervening SSRs. Variation in number of tandem 
repeats results in PCR product variation. The high degree of polymorphisms of SSRs results 
from high copy number of the basic motifs or internal heterogeneity of the sequences. 
Microsatellite markers has also been referred to as Simple Sequence Length Polymorphism 
(SSLP, Tautz 1989) Sequence tagged Microsatellite (STMS, Beckmann & Soller, 1990) 

The procedure for STMS primer directed PCR is also very simple. The region 
containing the microsatellite is amplified is amplified by PCR using primers that flank the 
microsatellite. The numbers of repeats present in the microsatellite or that allele determines 
the size of the DNA amplified. If the number of repeats of one allele were different from the 
other, than two separate bands would show up on the gel. 

A wide variety of methods for construction of libraries enriched for microsatellite 
sequences have been reported including those based on vectorette PCR using anchored 
primers (Lench et al, 1996), hybridization of adapter ligated DNA to the primer bound filters 
(Kandpal et al, 1994; Edwards et al, 1996), extension of uracil containing DNA templates 
using repeat primers (Ostrander et al, 1992) and hybridization of PCR generated single 
stranded template to biotinylated oligos (Kijas et al, 1994). Most of these involved in 
construction of either genomic or cDNA libraries in the Plasmid or phage vectors. 

Majority of approaches for the preparation of microsatellite markers is based either on 
the computer analysis of DNA sequence databases or the screening of various kinds of 
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genomic DNA libraries. In order to improve efficiency, microsatellite-enriched libraries have 
been constructed using various methods. Over the years the consensual approach has been 
developed that is largely based on the adapter-ligation method described by Kandpal, 1994 
and Hamilton, 1999. This involves shearing or sonication of genomic DNA to 400 bp- 
lOOObp sizes followed by ligation of universal adaptors. Then various di, tri or tetra repeat 
oligos are end labeled with biotin and hybridized to adaptor ligated size fractionated DNA. 
The corresponding is then captured using the streptavidin coated magnetic particles and 
separated on magnetic rack. These captured fragments are then eluted, washed and ligated to 
vector of choice and transformed to E. coli preferably using blue white selection. The 
transformants are then picked up on a nylon membrane and screened using same oligos that 
were earlier used for capture. These positive clones were sequenced subjected to various 
software’s like SSRIIT and TRF for analysis. Subsequently primers are designed flanking the 
repeat and used for amplification. 

Microsatellites are used as DNA markers for a variety of purposes. They can be used 
to construct genetic maps for e.g. to investigate QTLs for various traits. They may also be 
used to investigate genetic diversity, both within and between species. Microsatellites are 
suitable markers to investigate intra-specific gene flow and demographic parameters such as 
population size fluctuation. They may also be suitable for phylogenetic studies. 

STMSs are easier to use than RFLPs owing to the smaller amount of DNA required, 
higher polymorphism and ability to automate assays. These markers can be easily exchanged 
between researchers because each locus is defined by primer sequences. STMS assays are 
more robust than RAPDs are more transferable than AFLPs. They are now replacing RFLPs 
in genetic mapping of crop plants. A combination of STMSs and AFLPs is used to produce 
detailed genetic maps. 

Unlike PCR with arbitrary primers, sequence-tagged sites (STS) are primers that are 
based on some degree of sequence knowledge. These unique, sequence-specific primers 
detect variation in allelic, genomic DNA. STS have a particular advantage over RAPDs in 
that they are codominant, that is, they can distinguish between homozygotes and 
heterozygotes. They also tend to be more reproducible, because they use longer primer 
sequences. However, they have the disadvantage of requiring some pre-existing knowledge 
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of the DNA sequence of the region, even if only for a small amount. The investment in effort 
and cost needed to develop the specific primer pairs for each locus is their primary drawback. 
As with RAPDs, using PCR produces a quick generation of data and requires little DNA. All 
STS methods use the same basic protocols as RAPDs (DNA extraction and PCR) and require 
the same equipment. 


2.9 Gene flow in relation with Vigna spp. 

Recent studies by Yamamoto et al., (2005) have indicated that spontaneous 
interspecies crossing commonly occurs among vascular plants, and therefore that horizontal 
gene flow from transgenic plants into wild relatives is unavoidable. Few surveys, however, 
have been conducted to determine an actual flow frequency for individual plant species. The 
results indicate that gene flow actually occurs among V. angularis complex plants, and that 
transferred genes might be stably maintained by the offspring. The present observation 
cannot be directly applied to transgenic plants due to the lack of information as to whether or 
not artificially introduced genes behave as do native genes after horizontal movement. 
However, as gene flow is found to be inevitable in this species, specific assessment of 
whether or not the gene to be introduced confers a significant selective advantage to the host 
is critical for utilization of genetically modified plants in future. 



MATERIAL &METHODS 


3. MATERIALS AND METHODS 


Materials 

3.1 Plant materials 

To attain the objectives of mapping MYMV resistance with DNA markers, it was proposed 
to develop F 2 population from a cross between 1PU 982, an improved high yielding 
recommended variety but highly susceptible to MYMV infection (Figure 1), and Cuttack 
Local, a local landraces (LLRs) from the Cuttack district of Orissa known to be a low yielder 
but resistant to MYMV disease. The cross was done at Allahabad and the FI progenies were 
also planted there. The F2 generation was grown in NBPGR, New Delhi along with the 
parents. 

For diversity analysis a total of 375 individuals of 25 accessions of (ten from each of 
vigna mungo, vigna sylvestris and five from vigna hainiana ) were used. The sample details 
are in table 3.1. 


Table 3.1: Accessions of diversity analysis. 


Pop Id 

Species 

Acc. 

1 

Mungo Cv 

BBD-08-01 B 

2 

Mungo Cv 

BBD-14-01B 

3 

Mungo Cv 

BB-03-09 

4 

V.Hyniana 

BB-07-01A 

5 

V.Hyniana 

BBD-15-01A 

6 

V. Mungo 

BB 03-10 

7 

v.Siivistris 

BBL-63-2K 

8 

v.Silvistris 

BBL-76-2000 
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V.Silvistris 


V.Mungo 

v.Silvestris 


V.Silvestris 

V.Mungo 

V.Mungo 


V.Mungo 


V.Mungo 


V.Mungo 


V.Silvestris 


19 

V.Hainiana 

20 

V.Silvestris 

21 

V.Hyniana 

22 

V.Silvestris 

23 

V.Silvestris 

24 

V.Silvestris 


BBL-80-2K 


IPU-99/10 


BBL-50-2K 


BBL-40-2K 


Vamban 2 


PDU-1 


TMV-1 


PDU-1016 


UG-414 


BBL-55-2K 

KP3-Rajgarh 


IC-251-407 


V.Hainiana 


BB-2641 


BB-2642 


BB-2649 
































Methods 


3.2 Molecular marker analysis 
3.2.1 DNA isolation 

Extraction of total genomic DNA was carried out using CTAB method given by Saghai- 
Maroof et al. (1984) with minor modifications of Porebski et al. (1997) and Khanuja et al. 
(1999). Solutions, reagents and instruments are described in Appendix I. 

Five gram of young and healthy etiolated seedling material per accession were ground 
to fine powder in pre-cooled mortar and pestle in liquid nitrogen. Powdered material was 
quickly transferred to centrifuge tubes containing 15 ml pre-heated CTAB extraction buffer 
and vortexed to disperse the tissue powder. The tubes were then incubated at 65° C for 1 
hour, with intermittent shaking. At the end of the incubation, tubes were cooled to room 
temperature and 15 ml chloroform: isoamyl alcohol (24:1) was added, following which 
mixing of contents by inversion was carried out for about 10 minutes. The samples were then 
centrifuged at 25° C for 1 5 minutes at 12,000 rpm. The aqueous layer was transferred to fresh 
centrifuge tubes and the DNA was precipitated by adding 0.7 volume of isopropanol. DNA 
was spooled out carefully and excess chemical was drained out. DNA was placed in a 1.5 ml 
microtube, washed two times with 70% ethanol and dried under vacuum. DNA pellet was 
dissolved in minimum volume of T: E (10:0.1) buffer. 

The dissolved DNA was transferred to 2ml Eppendorf tubes; 2 pi of RNase A (25 
mg/ml) was added and incubated for 1 hour at 37° C. Subsequently, DNA was treated with an 
equal volume of phenol and spun in a microcentrifuge at 12,000 rpm for 15 minutes at 25° C. 
This step was repeated with phenol: chloroform (1:1 mixture) and chloroform: isoamyl 
alcohol (24:1 mixture). RNA-free, purified DNA was precipitated by adding 0.1 volume 
sodium acetate and 2.5 volumes chilled absolute ethanol to the aqueous phase and then 
collected by spinning at 12,000 rpm for 15 minutes. The precipitate was washed twice with 
70 % ethanol, air-dried and dissolved in TE buffer. 
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DNA concentration was determined by using DyNA Quant Fluorometer. Prior to 
quantification, the fluorometer was calibrated with 2 pi (1 mg/ml) of standard calf thymus 
DNA in 2 ml of low range assay solution (Appendix I). For quantification, 2 pi of DNA 
sample was added to 2 ml assay solution, mixed thoroughly and measurement in ng/pl 
recorded. 

Part of each DNA sample was diluted with sterilized distilled water to yield a 
working concentration of 20 ng/pl. The diluted samples were stored at 4°C for immediate use 
while the original concentrates were kept for long term storage at -20°C. 

3.22TDNeasy® Plant Miniprep Kit (Qiagen) 

To the 100 pi of grounded tissue added 400 pi of lysis buffer (API) and 4 pi of RNase 
solution (lOOmg/ml) and vortexed vigorously. Mixture was incubated at 65° C for 10 min 
and mixed 2-3 times during incubation by inverting tubes. To the lysate added 130 pi of 
precipitation buffer and incubated mixture for 5min on ice. Lysate was applied to QIA 
shredder spin column (lilac) placed in a 2 ml collection tube and centrifuged for 2min at 
14000 ipm. Flow through was transferred to a fresh collection tube. To the cleared lysate 1.5 
volumes of binding buffer (AP3/E) was added and mixed by pippetting. 650 pi of this 
solution was applied to DNeasy mini spin column sitting in 2ml collection tube and tubes 
were centrifuged for 1 min at 8000 rpm. Flow through was discarded. To the DNeasy mini 
spin column 500 pi of washing buffer (AW) was added and tubes were centrifuged for 1 min 
at 8000 ipm. Washing step was repeated once again. The DNeasy mini spin column was 
transferred to a fresh 1.5 ml microfiige tube and 100 pi of elution buffer (AE) applied to the 
column. Columns were incubated at room temperature for 5 min and centrifuged for lmin at 
8000 rpm to elute. The DNA was obtained as flow through. The eluted DNA was stored at - 
20°C 
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3.3 Random Amplified Polymorphic DNA (RAPD) analysis 
3.3.1 DNA Amplification 

Sixty six deca-nucleotide primers from kits A, B, C, D, M and O of Operon 
technology (USA) was screened for suitability of RAPD amplification. The nucleotide 
sequences of these primers and characteristics are listed in Table 2. The reaction mixture 
contained 3mM MgCl 2 , 0.2 mM dNTPs, 1U Taq (Bangalore genei), 0.4 mM primer and 20 
ng genomic DNA. The gene amp 9600 PCR system (Perkin Elmer Cetus, Norwak, CT) was 
used for amplification with following temperature profile: 94°C for 3 min; 40 cycles of 94°C 
for 1 min, 32°C for 1 min and 72°C for 1 min; 72°C for 5 min; 4°C end. 


Table 3.2: List of-sereened RAPD primers for mapping population 


S.No. 

Primer Name 

Primer Code 

Base sequence (S’ to 3’) 

1 

OPA2 

S22 

TGCCGAGCTG 

2 

OPA3 

S23 

AGTCAGCCAC 

3 

OPA4 

S24 

AATCGGGCTG 

4 

OPA5 

S25 

AGGGGTCTTG 

5 

OPA6 

S26 

GGTCCCTGAC 

6 

OPA7 

S27 

GAAACGGGTG 

7 

OPA8 

S28 

GTGACGTAGG 

8 

OPA10 

S30 

GTGATCGCAG 

9 

OPA11 

S31 

CAATCGCCGT 

10 

OPA13 

S33 

CAGCACCCAC 

11 

OPA14 

S34 

TCTGTGCTGG 

12 

OPA15 

S3 5 

TTCCGAACCC 

13 

OPA16 

S36 

AGCCAGCGAA 

14 

OPA17 

S37 

GACCGCTTGT 

IS 

OPA18 

S38 

AGGTGACCGT 
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16 

OPA19 

S39 

CAAACGTCGG 

17 

OPA20 

S40 

GTTGCGATCC 

18 

OPD20 

S60 

ACCCGGTCAC 

19 

OPD13 

S53 

GGGGTGACGA 

20 

OPB1 

SI 

GTTTCGCTCC 

21 

OPB2 

S2 

TGATCCCTGG 

22 

OPB3 

S3 

CATCCCCCTG 

23 

OPB4 

S4 

GGACTGGAGT 

24 

OPB5 

S5 

TGCGCCCTTC 

25 

OPB6 

S6 

TGCTCTGCCC 

26 

OPB8 

S8 

GTCCACACGG 

27 

OPB9 

S9 

TGGGGGACTC 

28 

OPB12 

S12 

CCTTGACGCA 

29 

OPB15 

S15 

GGAGGGTGTT 

30 

OPB16 

S16 

TTTGCCCGGA 

31 

OPB17 

S17 

AGGGAACGAG 

32 

OPB18 

S18 

CCACAGCAGT 

33 

OPB20 

S20 

GGACCCTTAC 

34 

OPC2 

S62 

GTGAGGCGTC 

35 

OPC4 

S64 

CCGCATCTAC 

36 

OPC6 

S66 

GAACGGACTC 

37 

OPCIO 

S70 

TGTCTGGGTG 

38 

OPC13 

S73 

AAGCCTCGTC 

39 

OPC15 

S75 

GACGGATCAG 

40 

OPC16 

S76 

CACACTCCAG 

42 

OPS0O 

SS01 

fiCTGGGGGAC 




43 

OPM2 

S402 

ACAACGCCTC 

44 

OPM3 

S403 

GGGGGATGAG 

45 

OPM4 

S404 

GGCGGTTGTC 

46 

OPM5 

S405 

GGGAACGTGT 

47 

OPM6 

S406 

CTGGGCAACT 

48 

OPM7 

S407 

CCGTGACTCA 

49 

OPM8 

S408 

TCTGTTCCCC 

50 

OPM9 

S409 

GTCTTGCGGA 

51 

OPM11 

S41 1 

GTCCACTGTG 

52 

OPM12 

S412 

GGGACGTTGG 

53 

OPM14 

S414 

AGGGTCGTTC 

54 

OPM15 

S415 

GACCTACCAC 

55 

OPM16 

S416 

GTAACCAGCC 

56 

OPM17 

S417 

TCAGTCCGGG 

57 

OPM18 

S418 

CACCATCCGT 

58 

OPM19 

S419 

CCTTCAGGCA 

59 

OPOI 

S44I 

GGCACGTAAG 

60 

0P02 

S442 

ACGTAGCGTC 

61 

OP03 

S443 

CTGTTGCTAC 

62 

0P04 

S444 

AAGTCCGCTC 

63 

0P05 

S445 

CCCAGTCACT 

64 

0P06 

S446 

CCACGGGAAG 

65 

0P07 

S447 

CAGCACTGAC 

66 

0P08 

S448 

CCTCCAGTGT I 




3.3.2 Preparation of agarose gel 

The amplification products were separated by electrophoresis in 1 .6% agarose gels according 
to the procedure outlined below. The gel trays were sealed by casting them in casting trays. 
The combs were placed in the slots and the whole set up was leveled. A 1.6% mixture of 
agarose was prepared by dissolving required quantity of agarose (3.84g) in 240ml of IX TAE 
buffer. It was boiled in a microwave oven till a transparent mixture was obtained. The 
mixture was cooled to about 50° C before pouring in gel trays. The trays were allowed to 
cool till the gel was firmly set. Then, placed in refrigerator at 4° C for 30 mins. 

3.3.3 Preparation of samples for loading in PCR tubes 

To each 25 ml mixture in PCR tubes, 2.5 ml of bromophenol blue-dye (10X) was added and 
the tubes spun for 5-10 sec. 

3.4 Sequence Tag Microsatellites Site (STMS) analysis 

For molecular characterization and population diversity analysis, 27 microsatellite primer 
pairs were used for amplification and resolution of alleles in the 25 vigna population. 

3.4.1 PCR amplification 

3.4. 1.1 Optimization of amplification conditions 

To determine optimal amplification reaction conditions, PCR was carried out at different 
concentrations of MgCk, Taq DNA polymerase (0.3 U and 0.5 U per 25 pi reaction mixture) 
and DNA (10 - 40 mg per 25 ml reaction mixture). The annealing temperatures tested were 
3°C above and below the melting temperature, T m of the particular primer. Conditions were 
finally optimized as given in Table 3.3, under which clear and consistent bands were 
amplified. A break-up of the primers based on repeat length and motif is depicted in Fig. 4. 
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Table 3.3: Primer details for STMS analysis. 


SL 

No. 

Primer 

Primer sequence (S’- 3’) 

Repeat Unit 

Annealing 
temp. (°C) 

1 . 

AB128093 

CCCGATGAACGCTAATGC TG 

(AG) 

57 



CGCCAAAGGAAACGCAGAAC 



2 

MB 122 A 

TGGTT GGTT GGTT C AC A AG A 

(TGGT) 

48 



CACGGGTTCTGTCTCCAATA 



3 

AB1281 13 

TCAGCAATCACTCATGTGGG 

(AG) 

55 



T GGG ACA AACCTC AT GGTTG 



4. 

AB128135 

ACTATTTCCAACCTGCTGGG 

(AG) 46 

48 



AGG ATTGT GGTT GGTGC ATG 



5. 

MB91 

GAGGCCAATCCCATAACTTT 

(AG) (GA) 

48 



AGC ACC AC AT C AG AG ATT CC 



6. 

VM22 

GCGGGTAGTGTATACAATTTG 

(AG) 

48 



GTACTGTTCCATGGAAGATCT 



7. 

VM24 

TCAACAACACCTAGGAGCCAA 

(AG) 

55 



ATCGTGACCTAGTGCCCACC 



8 

VM27 

GTCC A A AGC A A AT GAB GT C A A 

(AAT)(TC) 

(AC) ’ 

48 



T G A AT G AC AATGAGGGT G C 



9 . 

VM 31 

GT GTT CT AGAGGGT GTG AT GGT A 

(CT„ 6 

56 



CGCT CTTCGTT GAT GGTT AT G 



10 

M323b 

T GCTTCCTTTT GTCTG A GTT AG A A 

(GA) 

57 



T GACGGAGAGAG AGAGAG AGAG 






3.4.J.2 Thermal cycler conditions 


The PCR amplification was carried out in a DNA thermal cycler (BIOER ™), which was 
programmed as follows: Initial denaturation for 4 minutes at 94° C followed by 32 cycles of: 

1) Denaturation at 94°C for 30 sec, (2) Primer annealing at optimized T a of the particular 
primer (Table 3.3) for 1 minute, (3) Primer extension at 72°C for 5 minutes. A 10 min final 
extension at 72°C followed at the end of 32 cycles. 

3.4.2 Gel electrophoresis 

Gels of various concentrations ie. 5%, 6%, 7%, 7.5% and 8% concentration were prepared. 
For casting, the gel stand was assembled. Gel cassettes were prepared by inserting spacers on 
two parallel sides of the gel plates. The cassettes were inserted into the casting stand and 
screws were tightened. Polymerization mixture containing appropriate concentration of 
acrylamide : bisacrylamide solution, distilled water, ammonium persulphate, TEMED and 
resolving buffer was poured into the cassette with a syringe. Combs were placed on the open 
notched side of the gel cassette for forming wells. After polymerization, combs were 
removed and gels were pre run for 1 5 min at 100 V. 

3.4.3 Data analysis 

The amplification products were scored across the lanes comparing their respective 
molecular weights. Each band was treated as one STMS allele. Scoring of bands was done 
from photographs. Homology of bands was based on distance of migration in the gel. 

3.4.4 Population genetic parameters 

Out of the 255 samples, a minimal set of 225 were chosen for population genetics studies. 
These 225 mango landraces were representative of 6 contiguous geographic regions covering 
20 districts of 12 states and union territories of India (Table..). Observed and expected 
heterozygosities (Ho and He) were used to measure genetic diversity using POPGENE 
(available online at http://www.ualberta.ca/~fveh/) . Deficits in heterozygosity were measured 
using fixation indices, reduction in heterozygosity from the Hardy-Weinberg equilibrium 
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were tested. The test for genome-wide departures from neutrality, such as those produced by 
demographic processes, the Ewens-Watterson test of neutrality (Watterson 1978, 1986), 
with probabilities calculated based on both homozygosity and Fisher’s exact tests (Ewens— 
Watterson-Slatkin’s exact test; Slatkin 1994, 1996), was performed using the program 
Arlequin Ver. 2000 (Schneider et al. 2000). The Ewens-Watterson test enables the detection 
of deviations from a neutral-equilibrium model as either a deficit or an excess of 
homozygosity relative to the number of alleles found at a locus. 

3.4.5. Analysis of molecular variance (AMOVA) 

The analysis of the distribution of genetic variation across the germplasm groups established 
a priori was done using the AMOVA option of Arlequin software package 
(Arlequinsoftware, Schneider et al. 2000; available online a thttp://lgb.unige.ch/arlequin/) . 
Fixation statistics were produced for individual SSRs and groups of germplasm. The 
significance of the estimates was obtained through permutation tests, using 1,000 
permutations. 

3.5 STS Analysis 

Six sequence tag microsatellites loci, were used for STS analysis in selected vigna 
population. Primers were synthesized in order to amplify the target loci. PCR amplification 
was carried out in 25 pi reactions using 40 ng template DNA, 0.2pm each of the forward and 
reverse primer, 0.3 Units of Taq DNA polymerase and 0.2 mM dNTPs. The PCR 
amplification was carried out in a DNA thermal cycler (BIOER ™), which was programmed 
as follows: Initial denaturation for 15 minutes at 94° C followed by 25 cycles of: 1) 
Denaturation at 94°C for 30 sec, (2) Primer annealing at optimized T a of the particular primer 
for 1 minute, (3) Primer extension at 72°C for 5 minutes. A 10 min final extension at 72°C 
followed at the end of 25 cycles. Amplified products were initially checked by loading 2 pi 
PCR product on a 1 .5 % agarose gel. The amplified products were purified and sequenced by 
capillary sequencing. The sequence data was trimmed and edited using Edit with Ridom 
Trace Edit 1.1.0 Version, Ridom GmBH, Germany. Alignment of sequences in different 
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cultivars was earned out using ClustalX and finally neighbour-joining trees were generated 
based on the sequences. 


Table 3.4: Forward and reverse primer sequence details for STS 


SI. 

No 

Primer 

Primer Sequence 

(5’- 3’) 

STS Size 

(bp) 

Annealin 

g 

Temp. 

1 

BV1 65289 

ATTCATGATAACTCGTCGGATCGCA 

349 

55 



TCCAACTACGAGCTTTTTAACTGCA 



2 

OLM1 

ggtctgcgagctgttttt ggagaag 

342 

II 



GCAATTCCCTCCTCAGCT AAAAGTG 



3 

UDPGD- 

B3 

TGGT G A AG ATTT GCT GC ATT GGTGC 

202 

55 



AAGAAGAGGTTCTTGCCACGGCA 



4 

SUSY- 8 

TCGC A AT G A ACCAC AC AG ATTT CA 

488 

55 



GTCCAACCTTGCCATGGTGAAGATA 



5 

SHMT- 1 

ACC AC A ACT CAC AAGT C ACTTC 

570 

55 



TT GCTG AG A ACCT GCT GTT GGT ATG 



6 

RNAR- 8 

GTTT GGCAG ATT GTT GGGGT GA AG A 

417 

55 
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3.5.1 Genei Quick PCR purification kit 


To the PCR product was added 5 volumes of binding buffer and this mixture was loaded onto 
a spin column placed in afresh 2.0 ml collection tube. Column was centrifuged at lOOOOrpm 
for 1 min and flow through was discarded. Column was placed in a fresh collection tube and 
to it was applied 500pl of wash buffer I. Column was centrifuged at lOOOOrpm for 1 min and 
flow through was discarded. 700 pi of wash buffer II was loaded into column and was 
centrifuged at lOOOOrpm for 1 min. Flow through was discarded and column in collection 
tube was spun for 3-5min to remove any traces of buffer. Column was placed in fresh 2.0ml 
collection tube and to it was added 50pl of the elution buffer. The column was incubated at 
room temperature for 1 min and then centrifuged for lmin at 10000 rpm to elute. The DNA 
was obtained as flow through. The eluted DNA was stored at -20°C. 

3.5.2 Sequence Tagged Sites (STS) Analysis 

For the STS analysis 25 samples from each populations of Vigna species were taken. Ten 
primer pairs was screened and six gave good amplification. The product was purified either 
by using gel purification system or by PCR purification kit. The number of segregating 
sites(S), observed nucleotide diversity per site between any two sequences assuming that the 
sample is random (p), number of haplotypes (H), haplotype diversity(Hd), average number of 
pairwise nucleotide differences within population(K), were calculated using DnaSP version 
5. The genetic differentiation among the populations was calculated in terms of fixation 
index (Fst) that estimates diversity within a subpopulation with respect to total genetic 
diversity. In addition, average number of pairwise nucleotide differences (Kxy), nucleotide 
substitution per site (Dxy), and net nucleotide substitution per sit (Da) between populations 
were also calculated. The sequences were edited and the translated amino acids were aligned. 
The above parameters were also estimated on DnaSP. Phylogenetic analysis was performed 
by neighbour-joining (NJ) method with kimura 2-parameter distance matrix (Kimura, 1980) 
MEGA version 3.0 (Kumar et al 2004). 
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4. RESULTS 


4.1 RAPD analysis 

4.1.1 Morphological characterization/evaluation 

For mapping of MYMV resistance with DNA markers, a F 2 population grown in NBPGR 
fields, along with the parents from a cross between IPU 982, an improved high yielding 
recommended variety but highly susceptible to MYMV infection, and Cuttack Local, a local 
landraces (LLRs) from the Cuttack district of Orissa known to be a low yielding but resistant 
to MYMV disease. Individuals F 2 plants will be scored following an interaction phenotype 
(IP) scale of 0-9, where 0 and 1 stands for immunity and resistance, whereas 3, 5, 7 or 9 
denote various degrees of susceptibility (Table 4.1). 

4.1.2 Optimization of RAPD reaction 

Template concentration of 20 ng was found to be the ideal concentration as it gave maximum 
number of reproducible bands. Similar variation in the concentration of Taq DNA polymerase 
and MgCl 2 showed best results with 1U of Taq DNA polymerase and 3mM of MgCl 2 . 

4.13 Primer selection and survey 

Out of the 66 primers surveyed, 24 were selected for the analyses as others gave either 
suboptimal, indistinct or monomorphic amplification products (Figure 2). 





Figure 2: RAPD profiles of parents and F 2 populations, used for screening as 


indicated in table 3.2. 






















































































Legend (Figure 3) 

Figure 3: STMS profile of 25 accession of Vigw indicating the polymorphism 
existing for the STMS markers, generated with the primer pairs 
AB128135. The numbers Indicated in each lane correspond to the 
sequence of accessions as indicated in table 3.1, The latte marked Ml and 
M2 is the DNA molecular weight standard lOfhhase pair ladder of MBI. 
Fermentas, (USA) 





4.2 STMS ANALYSIS 


4.2.1 Primer screening 


A total of 30 STMS primers were screened to identify ideal primers for detailed population 
genetic studies of 25 populations of Vigna. DNA extracted from 15 individual plants per 
population was used for STMS analysis. The properties of STMS primers, primer sequence, 
repeat motif, their chromosomal location and annealing temperature are present in Table 2. 
Of the thirty primers screened, ten polymorphic primers were used in the present study for 
final STMS analysis. Annealing temperature ranged from 48°C to 55°C for the screened 
primers. 

Table 4.2 The number of alleles per locus and their allele size range 

Locus 

No. of 
alleles 

Allele size (bp) 

Range 

Difference 

AB 128093 

2 

180-200 

20 

MB 122 A 

3 

180-250 

70 

AB1281 13 

4 

150-170 

20 

AB128135 

3 

150-210 

60 

MB 91 

2 

150-180 

30 

VM 22 

3 

210-230 

20 

VM 24 

2 

150-170 

20 

VM 27 

5 

180-230 

50 

VM 31 

2 

160-170 

10 

MB323b 

4 

300-320 

20 

Total 

30 

*» 

— 
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Mean 


3.0 


4.2.2 Properties of STMS primers 

A total of 30 alleles were detected in 25 Vigna populations by using ten STMS markers. The 
number of alleles per locus and their size range are presented in Table 4. The number of 
alleles per locus ranged from 2 (AB128093, MB 91, VM 24 and VM 31) to 5(VM 27), with 
an average of 3.0 alleles per locus. The overall size of amplified products ranged from 150 bp 
to 320 bp. The representative amplification profiles are shown in Figure 3, 4 and 5. The size 
difference between smallest and largest alleles for a given STMS locus varied from 10 bp 
(VM 22 and VM 31) to 70 bp (MB122A). 

4.2.3 Genetic structure of Vigna populations 

The genotypes of fifteen individuals of each Vigna population were not identical at ten STMS 
loci which were studied. The allele frequency per locus generated by ten STMS primers in 25 
Vigna populations was 0.03 to 1.00 presented in Annexure II. The mean number of alleles 
per locus was 3.0 and that for group I ( Vigna mungo ), II ( Vigna mungo war. silvestris ) and III 
( Vigna hainiana) were 1.70, 2.90 and 2.20, respectively. The number of alleles per locus 
ranged from 2 (AB 128093, VM 22, VM 24 and VM 31) to 5(VM 27). Overall, a total of 30 
alleles were detected of which 8 were common (occurring with =0.05 frequency) and 22 were 
rare (occurring with <0.05 frequency). Of the 25 Vigna populations, for all the ten loci, the 
maximum number of 5 alleles was recorded for the populations BBL-50-2K (V. silvestris ) and 
minimum was two alleles for BB-07-01A ( V. hainiana). The number of alleles per locus 
ranged from 2 (AB128093, MB 91, VM 24 and VM 31) to 5(VM 27). Population wise the 
number of observed alleles per locus ranged from 1.0 (hainiana BB-07-01A) to 1.9 (V. 
sivestris BB-2642), The genetic diversity in terms of Shannon’s information index was 
highest (0.466) for the population BB-2642 ( sivestris ) and lowest (0.0) for the population BB- 
07-01 A (hainiana). The observed heterozygosity was highest (2.0) for the population 
IC251377 (hainiana) and lowest (0.0) for the populations BB-07-01A (hainiana). The 
expected heterozygosity was highest (0.32) for the population BB-2642 (sivestris) and a 
lowest (0.0) for BB-07-01 A (hainiana). The overall observed number of alleles and effective 
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Legend (Figure 5) 

Figure 5: STMS profile of 25 accession of I'igna indicating the polymorphism 
existing for the STMS markers, generated! with the primer pain MB 
122A. The numbers indicated in each lane correspond to the sequence of 
accessions as indicated in table 3.1. The lane marked Ml and M2 is the 
DNA molecular weight standard 100-base pair ladder of MB!. 
Fermentas, (USA) 




number of alleles were 3.0 and 1.7 with the genetic diversity in terms of Shannon’s 
information index was 0.62. The overall observed heterozygosity and expected 
heterozygosity were 0.08 and 0.35 respectively. The summary diversity analysis of the 25 
populations is presented in the Table 4.3. 

4.2.4 Genetic structure of grouped populations 

Groupwise allele frequency is presented in Table (Annexure II). Groupwise population 
genetic structure is presented in Table 4.3, 4.4 and 4.5. 

Group I (Vigna mango) contain ten populations (BBD-08-01B, BBD-14-01B, BB-03- 
09, BB_03-10,IPU-99/10,Vamban 2,PDU-1,TMV-1,PDU-1016 and UG-414). The number of 
alleles detected per locus varied from 1 (AB1281 13, MB 91, VM 22, VM 24 and VM 31) to 
4 (VM27). 

The mean number of observed and effective alleles was detected 1.7 and 1.1. The 
Shannon’s information index was 0.18 for the group 1. The observed heterozygosity and the 
expected heterozygosity was found 0.083 and 0.10 respectively. The F-statistics and gene 
flow for the group 1 was found 0.35 and 0.46. The genetic identity and distinctness among 
the populations of the group is given in the Table 4.6. Dendogram for the group 1 is given 
Figure 6. 

Group II (silvestris) contain ten populations (BBL-63-2K, BBL-76-2000, BBL-80-2K,BBL- 
50-2K,BBL-40-2K,BBL-55-2K,IC-25 1 -407,BB-2638,BB-2641 and BB-2642). The number 
of alleles detected per locus varied from 1 (VM 24) to 5 (VM27). The mean number of 
observed and effective alleles was detected 2.9 and 1.6. The Shannon’s information index 
was 0.57 for the group 2. The observed heterozygosity and the expected heterozygosity was 
found 0.094 and 0.34 respectively. The F-statistics and gene flow for the group 2 was found 
0.51 and 0.24. The genetic identity and distinctness among the populations of the group is 
given in the Table 4.7. Dendogram for the group 2 is given Figure 7. 

Group 3 ( Vigna huiniana) contain five populations (BB-07-01A, BBD-15-01A,KP3- 
Rajgarh,lC-25 1-377 and BB-2649). The number of alleles detected per locus varied from 1 
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Figure 4 : STMS prolife „r « 

“ iSti " 8 for «>' STMS n TrZl generated ,,,C P0 '™“ rphis 

The " Umbe « indicated in each 1 "" h Primer P a '> a VM2' 

ureessions as indicated in table 3" 1 7 » ** «,„en« . 

“ NA m ° lKUl - «i gh , slandard ^ "***"' and M2 is th. 
Fcrmentas, (USA) °0-iwse pair ladder of MBI 




(AB128093, VM 22 and VM 31) to 4 (VM27). The mean number of observed and effective 
alleles was detected 2.2 and 1 .7. The Shannon’s information index was 0.52 for the group 3. 
The observed heterozygosity and the expected heterozygosity was found 0.06 and 0.31 
respectively. The F-statistics and gene flow for the group 3 was found 0.76 and 0.07. The 
genetic identity and distinctness among the populations of the group is given in the Table 4.8. 
Dendogram for the group 3 is given Figure 8. 

Overall in three groups the allele frequency was varied from 0.01 to 1.00 and found high in 
group 2 than other two groups but observed number off alleles and effective number of 
alleles was high in group2 (2.9) and group3 (1.7). The observed heterozygosity (0.09) and 
expected heterozygosity (0.33) was high in group2. The F-statistics and gene flow was high 
in group3 (0.76) and group 1 (0.4647) respectively. 
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Table 4, 3 Summary diversity of 25 Vigna populations based on STMS primers 
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Table 4A Administrative groupwise summary genetic variation statistics for all loci of 25 Vigna populations 
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Table 4.6 Nei’s Unbiased Measures of Genetic Identity and Genetie distance (Group 1). 


12 3 6 


1 

sjesieslcsic 

0.980 

0.975 

0.952 

2 

0.020 

**** 

0.977 

0.931 

3 

0.026 

0.023 

**** 

0.927 

6 

0.049 

0.071 

0.075 

**** 

10 

- 0.000 

0.022 

0.026 

0.048 

13 

0.005 

0.017 

0.035 

0.056 

14 

- 0.000 

0.015 

0.027 

0.050 

15 

- 0.000 

0.022 

0.026 

0.048 

16 

0.105 

0.044 

0.138 

0.161 

17 

0.065 

0.022 

0.097 

0.120 


10 

1.000 

0.978 

0.974 

0.953 

**** 

0.008 

0.001 

- 0.001 

0.110 

0.070 


13 

0.995 

0.983 

0.966 

0.945 

0.992 

0.003 

0.008 

0.084 

0.051 


14 

1.000 

0.985 

0.973 

0.952 

0.999 

0.997 

**** 

0.000 

0.090 

0.054 


15 

1.000 

0.978 

0.974 

0.953 

1.001 

0.992 

0.999 

0.110 

0.070 


16 

0.901 

0.957 

0.871 

0.852 

0.895 

0.920 

0.913 

0.895 

**** 

0.030 


17 

0.937 

0.978 

0.907 

0.887 

0.933 

0.951 

0.947 

0.933 

0.997 

**:*:*: 


Figure 6. Dendrogram based on Nei’s (1972) Genetic distance: Method = UPGMA 

Modified from NEIGHBOR procedure of PHYLIP Version 3.5 (Group I) 
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Table 4.7 Nei's Unbiased Measures of Genetic Identity and Genetic distance (Group 


1 

2 

3 

6 

10 

13 

14 

15 

16 
17 


1 2 


**** 

0.980 

0.020 

**** 

0.026 

0.023 

0.049 

0.071 

- 0.000 

0.022 

0.005 

0.017 

- 0.000 

0.015 

- 0.000 

0.022 

0.105 

0.044 

0.065 

0.022 


3 6 

0.975 

0.952 

0.977 

0.931 

**** 

0.927 

0.075 

**** 

0.026 

0.048 

0.035 

0.056 

0.027 

0.050 

0.026 

0.048 

0.138 

0.161 

0.097 

0.120 


10 13 

1.000 0 . 9 ^ 

0.978 0.983 

0.974 0.966 

0.953 0.945 

**** 0.992 

0.008 **** 

0.001 0.003 

- 0.001 0.008 

0.110 0.084 

0.070 0.051 


14 ' 

15 

1.000 

1.000 

0.985 

0.978 

0.973 

0.974 

0.952 

0.953 

0.999 

1.001 

0.997 

0.992 

**** 

0.999 

0.000 

**** 

0.090 

0.110 

0.054 

0.070 


16 

17 

0.901 

0.937 

0.957 

0.978 

0.871 

0.907 

0.852 

0.887 

0.895 

0.933 

0.920 

0.951 

0.913 

0.947 

0.895 

0.933 

**** 

0.997 

0.030 

**** 


Figure 7 Dendrogram based on Nei’s (1972) Genetic distance: Method = UPGMA 
■ Modified from NEIGHBOR procedure of PHYLIP Version 3.5 (GroupH) 
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Table 4.8 Nei’s Unbiased Measures of Genetic Identity and Genetic distance (Group 


4 

5 

19 

21 

25 


0.9377 

0.7483 

0.4748 

0.6179 

0.0643 


0.7180 

0.4984 

0.6837 

0.2900 

0.3313 


0.7672 

0.7561 

0.7449 

0.6963 

0.2649 


0.5396 

0.4815 

0.3802 

0.2796 

0.6169 



Figure 8 Dendrogram based on Nei’s (1972) Genetic distance: Method = UPGMA 
Modified from NEIGHBOR procedure of PHYLIP Version 3.5 (Group III) 
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4.5 Evidence of group differentiation and selection 

Group differentiation is indicated by F ST value (Table 4.5), which measures the 
fixation of alleles in different populations. F ST values ranged from 0.35 to 0.76 this 
indicates that there is moderate to high genetic differentiation between populations. 

Ewens-Watterson test for neutrality was done using 1000 simulated samples 
(Table 4.9). All the loci were neutral to selection pressure as their observed F value was 
between lower and upper 95% limit. 

Table 4.9 The Overall Ewens-Watterson Test for Neutrality for all loci. 


s. 

No. 

Locus 

Obs F 

L 95 

U95 

1 

AB 128093 

0.84 

0.51 

1.00 

2 

MB 122 A 

0.44 

0.39 

0.99 

3 

AB1281 13 

0.54 

0.35 

0.98 

4 

AB128135 

0.71 

0.39 

0.99 

5 

MB91 

0.56 

0.51 

1.00 

6 

VM22 

0.95 

0.40 

0.99 

7 

VM24 

0.73 

0.50 

1.00 

8 

VM27 

0.41 

0.30 

0.97 

9 

VM31 

0.92 

0.50 

1.00 

10 

MB323b 

0.34 

0.34 

0.98 


Major genetic differentiation was found between wild and cultivated species of 
Vigna. Using POPGENE software, grouped the populations into two broad clusters. The 
first cluster (Population BBD-08-01B, BBD-14-01B, BB-03-09, BBJB-10, IPU-99/10, 
Vamban 2, PDU-1, TMV-1, PDU-1016 and UG-414) comprising cultivated varieties of 
Vigna mungo as well as 10 population of (BBL-63-2K, BBL-76-2000, BBL-80-2K, 
BBL-50-2K, BBL-40-2K, BBL-55-2K, IC-251-407, BB-2638, BB-2641, BB-2642) wild 
relative Vigna sylvestris. With in this cluster the cultivar formed an independent 
subcluster, closer to subcluster II consist of six populations of Vigna sylvestris. Four 
population of Vigna sylvestris (20, 22, 23 and 24) formed subcluster III. Wild population 
of Vigna hainiana (five populations) formed a distant cluster (Figure 9).The genetic 
identity and distinctness among the populations of the group is given in the Table 4.10. 
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able 4.10 Nei's Unbiased Measures of Genetic distance for all the Vigna populations 
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Figure 9 Dendrogram based on Nei’s (1972) Genetic distance: Method = UPGMA 
-Modified from NEIGHBOR procedure of PHYLIP Version 3.5 
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Figure 10: STS profile of 25 accession of Vigna indicating the screening of STS 
markers. The lane marked Ml and M2 is the DNA molecular weight 
standard 100-base pair ladder of MBI. Fermentas, (USA) 







4.6 Analysis of molecular variance 


For realizing the partitioning of molecular variance at different levels of hierarchy. 
Analysis of molecular variance (AMOVA) was performed using Arlequin software. 
AMOVA revealed maximum percentage of variation within populations (16 03%) 
followed by 39.13% among populations within group and 44.84% among groups (Table 
4.1 1). Average F-Statistics over all loci were, FSC: 0.710, FST: 0.840 and FCT: -0.840. 

Table 4.11 Analysis of molecular variance (AMOVA) at different geographic 
levels based on Arlequin 


Source of 
variation 

Sum of 
squares 

Variance 

components 

Percentage 

variation 

Among 

groups 

333.49 

1.25 

44.84 

Among 

populations 

within 

groups 

369.67 

1.10 

39.13 

Within 

populations 

156.40 

0.45 

16.03 

Total 

859.56 

2.79 



4.7 Sequence Tag Microsatellites Analysis 

The present study was undertaken to assess genetic diversity and relationships among 
selected accessions. Twenty five samples were taken from each population. A total of 
twelve primers was screened and out of twelve six primers gave reproducible bands and 
were studied for analysis. Amplified products were purified using PCR product 
purification kit (Banglore Genei). Purified products were checked by running 2 pi of 
amplified products on 0.8% agarose gel (Figure 10). 

4.8 Sequencing 

Purified products were sequenced directly using automated DNA sequencer, ABI GA310. 
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4.9 DNA sequence analysis 

The DNA sequence data was first edited using Ridom Trace Edit version 1.1.0 Ridom 
GmbH, Germany. The clarity of peak separation and height were the criteria used for 
ascertaining the quality. The acceptable sequence traces were included for further 
analysis. The sequences were then aligned using CLUSTAL X software (Thompson et 
al., 1 997). The gap penalties and weights were adjusted to obtain reliable alignments. The 
aligned sequence matrices were used for phylogenetic tree reconstmction using the 
maximum parsimony, minimum evolution and neighbor joining criteria with MEGA 4.1 
software (Kumar et al., 2008). The statistical support for the branches was determined by 
bootstrap analysis of 5000 psedo-data sets. The percent bootstrap values were indicated 
on each of the branch point on the trees. 


Table 4.12 Measures of DNA sequence polymorphisms and tests of neutrality 


Locus 

Segregating 

site 

Singleton 

Variable 

sites 

Parsimony 

Informative 

site 

Total no 

of 

mutation 

K 

H 

Hd± 

S.D 

p ±S.D 

Tajima’s 

D 

RNAR8 

16 

9 

1 

16 

2.58 

5 

0.967± 

0.019 

0.01 299± 
0.00481 

-1.639 

SHMT 1 

3 

3 

0 

3 

0.40 

3 

0.257 ± 
0.142 

0.00094 

± 

0.00058 

-1.685 

SUSY-8 

11 

3 

4 

7 

1.24 

4 

0.333 ± 
0.124 

0.00907 

± 

0.00339 

-1.130 

BV16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

UDPGDB3 

5 , 

2 

3 

5 

0.96 

3 

0.338 ± 
0.120 

0.00982 

± 

0.00365 

-0.928 

OLM1 

47 

5 

42 

53 

17.37 

14 

0.971 ± 
0.032 

0.07297 

± 

0.00478 

0.569 


S; Number of segregating (polymorphic/variable) sites; K; Average number of pairwise 
nucleotide differences, H; Number of haplotypes, Hd; Haplotype diversity; p; Observed 
average pairwise nucleotide diversity; D; Tajima’s D test statistics. 

The measures of DNA sequence polymorphism and tests of neutrality was found 
very high in OLM1 other than five loci. The number of haplotypes in locus OLM1 (H = 
14) was higher than RNAR8 (H= 5), SHMT 1 (H= 3) and SUSY 8 (H= 4). The number 
of segregating site, parsimony informative site and total number of mutations were 
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highest 47, 42 and 53 m locus OLM1. Singleton variable site lound highest y in locus 
RNAR8. The haplotypes (gene) diversity (Hd) was found highest 0.971± 0.032 SD. The 
Average number of pairwise nucleotide difference (K) was 17.37 with overall p diversity 
0.07297 ± 0.00478 SD in the locus OLM1 (Table 4.12). The positive value 0.569 for 
Tajima’s D indicates positive natural selection. 

Inter population nucleotide difference (Kxy) varied from 0.0235 (between 
silvestris and mungo) in locus UDPGDB3 to 26.639 (between mungo and hainiana ) in 
OLM1. Similarly, the average number of nucleotide (Dxy) and net nucleotide (Da) 
substitutions per site between populations ranged from 0.00088 (SHMT1) to 0.09146 
(OLM1) respectively. Inter-population comparison of the total (TM) versus shared (SM) 
number of mutation showed only on locus OLM1 (29 of 53 mutations were shared) 
(Table 4.13). Gene flow was varied from 0.31 (SUSY8) to 8.83 (OLM1) on basis of 
haplotype data information (Table 4.14). Based on sequence data information it was 
found 0.15 (SUSY8) to 7.1 1 (OLM1). 

Table 4.13 Inter-population genetic differentiation 


S. Population-] 
No 

1 mungo 

2 mungo 

3 silvestris 

4 mungo 


Population-2 Kxy 
silvestris 0.71 ' 


0.00369 


Da TM/ 

SM 

0.00022 4 


hainiana 


hainiana 


0.02740 0.01461 21 
0.03551 0.01865 14 


mungo 


silvestris 


SHMT 


7 hainiana 


8 hainiana 


mungo 


10 silvestris 


1 1 silvestris 


12 mungo 


1 3 mungo 


silvestris 17.375 0.07445 0.00309 53/29 

hainiana 26.639 0.09146 0.05232 72 

hainiana 15.028 0.06015 0.02832 34 

mungo 4.045 0.02173 0.01238 16 

silvestris 2.076 0.02676 0.01622 7 

silvestris 0.444 0.00324 0.00000 4 

mungo 0.235 0.00255 0.00000 2 

hainiana 1.561 0.02525 0.01515 5 

hainiana 1.154 0.02273 0.01515 3 


silvestris 0.400 0.00088 0.00085 3 


Kxy; Average number of nucleotide differences between populations 


% 


I 


Uxy; The average number ot nucleotide substitutions per site between populations. 

Da; The number of net nucleotide substitutions per site between populations, 

TM, total number of mutations in both populations; SM, Numbers of mutations shared 
between populations 


Table 4.14 Gene Flow Estimates based on sequence analysis.(Nei 1973) 



Haplotype Data 
Information 

Sequence Data Information 



Primer 

GST 

Nm 

DeltaSt 

GammaSt 

Nm 

Fst 

Nm 

SHMT1 

0.03493 

6.91 

0.00006 

0.06250 

3.75 

0.501 

0.25 


0.25064 

0.75 




0.00 



0.02752 

8.83 

0.00584 

0.08384 

2.73 

0.034 

7.11 

SUSY8 

0.44439 

0.31 

0.00514 

0.57108 

0.19 

0.628 

0.15 

~UDPGDB3 

0.42677 

0.34 

0.00546 

0.58416 

0.18 

0.600 

0.17 


Table 4.15a Population wise Measures of DNA sequence polymorphisms and tests of 
neutrality for RNAR8. 


S. no 

Population 

5 

h 

Hd 

K 

Pi 

JC, PiJC 

1 

mungo 

2 

2 

0.33333 

0.66667 

0.00335 

0.00337 

2 

silvestris 

2 

2 

0.35556 

0.71111 

0.00357 

0.00360 

3 

hainian 

9 

2 

0.66667 

6.00000 

0.03015 

0.03110 


Total 

16 

5 

0.53216 

2.58480 

0.01299 



Number of segregating sites, S: Number of haplotypes, h: Haplotype diversity, Hd: 
Average number of differences, K: Nucleotide diversity, Pi: Nucleotide diversity with JC, 
PiJC: 
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■ able 4.15b UDPGDB3 


S. no 

Population 

S 

h 

Hd 

K 

Pi 

JC, PiJC 

i 

silvestris 

2 

2 

0.25000 

0.50000 

0.00510 

0.00517 

2 

mungo 

0 

1 

0.00000 

0.00000 

0.00000 

0.00000 

3 

hainiana 

3 

2 

0.50000 

1.50000 

0.01531 

0.01563 


Total 

5 

3 

0.33810 

0.96190 

0.00982 



Table 4.15c SUSY8 


S. no 

Population 

S 

h 

Hd 

K 

Pi 

JC, PiJC 

1 

hainiana 

3 

2 

0.66667 

2.00000 

0.01460 

0.01482 

2 

mungo 

0 

1 

0.00000 

0.00000 

0.00000 

0.00000 

3 

silvestris 

4 

2 

0.22222 

0.88889 

0.00649 

0.00662 


Total 

7 

4 

0.34762 

1.29524 

0.00945 

— 


The total no of segregating sites (S) for Population wise Measures of DNA 
sequence polymorphisms and tests of neutrality was varied from 16 (RNAR 8) to 5 
(UDPGDB3). Number of haplotypes also high 5 in RNAR 8 where for the rest two locus 
were 3 (UDPGDB3) and 4 (SUSY8) respectively. Haplotypes diversity and nucleotide 
diversity value were 0.53216 and 0.01299 which is very high value than other two loci. 
Individually Vigna hainiana showed maximum value in all respect for all the three loci 
(Table 15a, 15b and 15c). 

4.10 Gene Tree 

The Minimum Evolution tree (Fig. 10a) for the locos BV 165289 had silvestris BB2638 a 
separate branch and silvestris BBL 63 2K and silvestris BBL 63 2K were much closed to 
each other. Neighbor Joining tree (Fig. 10b also gave the same pattern as ME tree only 
species were different. 
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According to CLUSTERX (1.83) MULTIPLE SQQUENCE ALIGNMENT in the locus 
BV1 65289 there was a deletion of a single nucleotide in silvestris BBL63-2K, BBL55- 
2K, BBL76-2K and mungo BBD 2649. There was an insertion of a single nucleotide for 
hainiana BB 2649. 

The minimum evolution tree (Fig. 11 a) for the locus OLM1 there was two main 
branches. In cluster I silvestris BBL76 2K was connected with mungo Vanban 2 with a 
bootstrap value of 93%. On another cluster, mungo BBD-08-01B, mungo BBD- 14-0 IB 
and mungo BB-03-09 were lying in a same group with bootstrap values of 100 %. In 


cluster II mungo IPU 99 was separate from all other with a bootstrap value of 30%. 
Neighbor Joining tree (Fig. 11 f) also gave the same pattern as ME tree. In Maximum 
parsimony (MP1) number of nucleotide insertion was three and the nucleotide was ‘A’ 
where as MP 2 number of insertion two (‘T’) in MP 3 and in MP 4 four (‘A’). For the 
locus SHMT1 the cluster was divide into two main clusters one was having only one 
species (silvestris BBL 40 2K) and another cluster containing rest all the accessions 
except silvestris 


Legend {Figure if) 


Figure II: Gene trees showing minimum evolution and 
locus BVI6 (#) ME, fb)NJ 


neighbour joining tree for 
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Figure 10(a) ME 
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Figure 10(b) NJ 
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Legend (Figure 12) 

Figure 12: Gene trees showing neighbour joining minimum evolution and 
maximum parsimony tree for locus OLM1 (a~f) 
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Figure 11(a) ME 


(b) MP1 


G silvestris BB2642 
G hainiana BB2649 
G silvestris BB2S38 
G silvestris BB2M1 
G mungoPU9910 
A mungo BBDU01B 
■A mungo BB03 09 
A mungo BBD0801B 
G silvestris BBLS02K 
G mungo PDU101S 
G silvestris BBL552K 
G mungo PDU1 
■G silvestris BBL802K 
G silvestris IC251407 
G silvestris BBL762K 
G mungo Vamban2 
G mungo UG414 




A silvestris BB2642 
■A hainiana BB2649 
•A silvestris BB2131 
•C silvestris BB2641 
■C mungo PU99 10 
A mungoBBDHOIE 
■A mungo BB03 09 
A mungoBBDOSOIE 
A slvestrisBBL50 2i 
A mungo PDU1018 
■A silvestris BBL55 21 
■A mungo PDU1 
■A silvestris BBUO 21 
•A silvestris K25UC7 
■A silvestris BBL7E 21 
-A mungo Vamban2 
•A mungo UG414 


(c) MP 2 


(d) MP 3 



G 


G 


(e) MP 4 


— t 
G 


G 


G 


G 




silvestris BB2642 
hainiana BB2649 
silvestris BB2638 
silvestris BB2641 
mungo IPU99 10 
mungo BBD14 01B 
mungo BB03 09 
mungo BBD08 01B 
silvestris BBL50 2K 
mungo PDU1016 
silvestris BBL55 2K 
mungo PDU1 
silvestris BBL80 2K 
silvestris IC251407 
silvestris BBL7 6 2K 
mungo Vambanl 
mungo UG414 






100 

24 


40 

- — - — silvestris B&2G41 


28 



silvestris BB2&38 


82 

, silvestris BB2642 




HmMmnm B&2649 


mungo BBD14 01 B 
mungo BBDOS 01B 
mungo BB03 09 

mungo 9*1199 10 


0 01 


(0 NJ 



IC251407. Neighbor Joining tree (Fig. 12b) also gave the same pattern as ME tree only 
species were different. 

Analysis of CLUSTERX in three species ( silvestris BB 264, BBL80-2K and mungo PDU 1) 
there is single nucleotide difference and SNP found in silvestris BB2641 and mungo IPU 
99/10. There was several nucleotide differences between accessions. 

Neighbor Joining tree and minimum evolution (ME) tree for the locus SHMT1 was mirror 
image of one another. There were two main cluster, one cluster consist only one accession 
silvestris BBL40-2K and another cluster consist rest of the samples. In CLUTERX deletion 
of single nucleotide found in silvestris BBL40-2K, BBL50-2K and mungo BB03-09 and 
insertion of single nucleotide in silvestris IC251407, BBL40-2K, BBL80-2K and mungo IPU 
99/10. Single nucleotide polymorphism (SNP) found in silvestris IC 251407 and BBL40-2K. 

Maximum evolution (ME) tree for the locus SUSY 8 had main two clusters. One consist only 
one accession silvestris IC251407 which is diverse from all others accessions. Cluster II 
divided into three sub cluster. Sub cluster one no groping found but in subcluster II hainiana 
BB0701A and hainiana BBD15-01A grouped together. Neighbor Joining tree (Fig. -) also 
gave the same pattern as ME tree only species were different. Analysis of CLUSTERX there 
was a single nucleotide deletion aid insertion found in hainiana BB07-01A and hainiana 
BB2649 respectively. Single nucleotide polymorphism (SNP) was detected in silvestris 
BB76-2O0O and IC 251407. 

Neighbor Joining tree and minimum evolution (ME) tree for the locus UDPGDB3 found only 
two cluster. Hainiana K.P3 Raj, IC25 1 377 and BB07-01 A were grouped together and rest 
all the accessions in another group. For the Maximum parsimony no such insertion was 
found. In CLUSTERX no deletion was found only insertion was found in silvestris BB2642. 
Single nucleotide polymorphism (SNP) was found only in hainiana BB2649. 
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Legend (Figure 13) 

Figure 13: Gene trees showing neighbour joining and minimum evolution tree for 
locus SHMT (a and b) 


Legend (Figure 14) 

Figure 14: Gene trees showing neighbour joining, minimum evolution and 
maximum parsimony tree for locus SUSY (a -c) 
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Legend (Figure 15) 

Figure 15: Gene trees showing neighbour joining minimum evolution and 
maximum parsimony tree for locus UDPGDB3 (a-c) 



silvestris BB2642 



mungo BBD14 01B 
mungo UC414 
mungo PDU1016 
silvestris BBL50 2K 
mungo PU99 10 
silvestris BB2S41 
mungo BBD08 01 B 
mungo PDU1 
mungo BB03 10 
silvestris BBL552K 
silvestris BBL632K 
mungo BB03 09 
silvestris $8178 2000 
silvestris BB2842 
mungo Vamban2 
silvestris BBL802K 
silvestris BBL40 2K 
mungo TMV1 
heinm&BQOJm 
twintem B8D15 01 A 

silvestris 1051407 


silvestris BB2638 
mungo TMV1 
mungo BB03 10 
mungo BBD14 01B 
silvestris BBL50 2K 
silvestris BBL55 2 K 
silvestris IC251407 
mungo BB03 09 
mungo BBD08 01B 
halniana BB2649 
mungo PDU1Q16 
silvestris BBL40 2K 
mungo PDU1 
mungo UG414 
silvestris BBL80 2K 
silvestris BB2641 
mungo Vamban2 


hainlana KP3 Rajgarh 
hainiana IC251377 
hainlana BB07 01 A 


0.002 


0.002 


(c) NJ (SUSY) 


Figure 14 (a) ME (UDPGDB3) 



s 

aitmtrisKZSim 

8 



mungo VG414 


mwgcPDtn 


sttvmis8SL402K 


sfmrn$&2m 


mvrnkmMm 


Mtveetrif B8L3Q 7K 


mtng&BBCWOlB 


mungo V&mbaff2 






wmgoBBtom 

j 

mmgo&mw 

1 

! 


| 

mmgoBmUMB 

i 

rnmg&mtm 

i 

mrngoWm 

1 



silvestris BB2W2 


h&mima KP3 Ra^gerts 

rrr 

**WR$*PI w # 




(b)M 



CLUSTAL X (1.83) MULTIPLE SEQUENCE ALIGNMENT 
m and Settings fCV.BhatMy DocumentsKoleKOLE SEQKDDPCDBffip.e: Thu Jun 04 10:2. « 2009 


Page 1 of 1 


flffffK 


Pg*;; mtSTt t!ih< 







m mpffi 

law^rt* 

ifaj&i mMni** 

M&fi^ 

Nffjif* MftW Wfef*» 

imit&i 

i&f** 

W^ffm 


mmi ii uu 

3* ? * ^ *wf? 
^ » ’i * *r*s 

« ? $ * WJ 

,ffl£ W 4» Wffff 

$ * i * *r*$ 
ft h% * 

ft 5 $ * >!Wff 
a** i *mr? 

~ a * * • < “ 

ft *** * «r v 

ft m *4 ^4 ftfff* 

ft w ft # *»«§#f* 

ft n '« * *"***; 

ft '4, » « pwv#* 

ft % £ * 1ft# 1 ** 
ft % t . * ftffft 
ft * • - •ft** 
mm .** m ** 

MM <*i M * f ‘" J "'‘ ‘ 1 

• » > 



|’"T 





lit 

m 

g n n 
in-n 
An n 
inn 

& 3? I 

in 


m •* 

f« *r 

' ]t;f r 


.: -; ,*«*>■ 
gf-ypa*. 


U'-/* 


Wffjfff' ftfgff Mjj*«? If Iff* if 4 ,*5fe J» 

$gt if n ^ fj ** * 

KWHtoftt ^ i*#*!* f < V. > it>y*t x m >* 

ewSfF. fff *■»?•!• >1 I f ! « r *> ** * 

jjpfipfff# ft* t - * * f5Y< ,l ^ !l ‘^? ? 

«rww v** w? ? * t i: ^ * r -w » 

fffftff *f ^ • “*v\ * j»* ?* 

yffttl" 6 ,1V r ; * ‘i h */,*v ^ «i? » 

^,/v Mf p % r * ti « in I 

y* f;i( »v^| i^r;j Ws t 4 4/ 4 * *»’ 1* *v l‘ 

^ff*fV ! ^ »« * »' ‘ ' ' 

fe|SisfS|f A 4? § % jmffl wf i* * itsft h 

frm <Hif^ #*■ 4 •< ' ' fV 1 ' ^ « 4 *t j. 

fpv \ * ‘ >' '* «» *.< 

Jg'flWV '*?•* J *•*»" j*» > 4 -< r i h di > 

*i‘«i-i ^wr* m f • \ l 1 0 '< i «*» 


• •:;r' l i S 

r , ,-l uM ■ i 

i | ©S I 

* » 

pm j 

v , J : 1 * 
x: f %,} f j 
U 1. - u 
4 « IflS 


» * i - >« 

f?’ "'il 

- ;' i ¥ ( > 
vi v » 5 ’> 

t ~ 

4^4. „(, 1 »* : 

'’*?'? ’ ; ?/ ; 
xpn\ rt t? i 


lrastri 8_BBL50_2K - 
mungo~IPU99__l 0 - 
lv^triflJB8ti7fi__23C - 
lvaBtris__BBL0O_2K 

mungo Vaab«n2 - 
lvastrlfl IC251407 - 

mungo PEUl - 
lvastriB__BBL55_23C - 
nmngo_D3414 - 
rular i 


mungo_BBDX4_01B 
mungo_BB03_09 
mungo BBD08 01B 
silvastr i a_BB?44 2 
hainian«~BB2 64 9 
BilveBtria“BB2«38 
flilvaetr i a~BB2€4 1 
mungo PDU1C15 
lvastri #_BBLS 0_2K 
mungo_IPU99_lO 
lvaatni#~8BI»76 2K 
lvastri b~BBL 80"2X 
mungo VesbSEn 2 
lvastri s_ > IC2514 07 
mungo PBU1 
lvastris_BBLS5^2K 
mungo_00414 











h'hgs® 

Lb.w..W k km' 


V ' i 

v* 

/ T fcfc 

J; A* .SfU ' * 

? ' ‘V® r 


' . " j P s 

& \ kf* r 
f & f ft 

»t [ M 

i » fft? 


f *¥& * i, > 4 -^1! * 

- mom* » * m 


« ***? % sg^f? r » tayi § 

", ff ,, ’Wm i> «. 5 * 

(tt Sf«jj!» -f- 3 i» «•% Mf. 

' « ,jy K « ^4r r <i «i»<m >'}Sk 

* * Hj|| I# A SW 1 ’ £ 

* > ; pi % >a ^ . - 5* a i j £» 

V ■". : :■"' :■ -.' ,: :! ■ . ■ 

* ' - ■ . »,.,•*-** jj* ’ . Mi« , k, : 

>- ! ■#,’1.' ^ it, 

f « V ’*<“$ - ' !» ",H 

; : 

' '■■" ■ ■■■■•. 

m .Kf^e ■'■? ■■* 1 ■; V; , ;:-z j:: . 


^ ■ fil '* ' r L/- , ' ; 

<Ws** «t*>rr(p <# «!':,, r«~t iti. ,» ,« *9 « 

: , ', , , 

I 

' . i ; 

... 

'- . :. 
i I - I f ', 1 ,, J, !„ u i 1 

:■•■;■■■ ■■■;■ 

4® sW*^ -/#■&?',,* |L i ;<« 2 A 

J i r !( "" t < ^ ’ *■ 

i) ; f 4 * JfisN; r , p 4 ’• 

' ■■ •; V' - ■ 

,■;■ . , 

■ ‘ ; 1 ’ 

j h f A *, ‘ If 


- !", ;;2 
: ; 


& 

>'* ; W 

i» fi. 


.ja 1 ^ , 2 , < 




Ivo&tris BBL&0J2K 

ivwtrielBBtH 2K 


I«.tri»jc25l407 

®ungo_M01 

lv«*txi*_BBI»55__2K 


jQungo BBD08J1B 
ailv*6trlfl_BB2642 
h4inian*_BB2649 
silv«etri6_BB2638 

gilY*»triB_BB2641 

jaungo_PDUlOlo 

lv«#tri«_BBL50_2K 
tnungo_IPU99 10 
lv«»tri*_BBIi7 6™2K 
lv«tri»_BBL80_2K 
Bungo_Vambaii2 
lv«trisJC2514t)7 
mungo_FDOl 
lv«stris_BBL55_2K 
nungo 00414 


CLUSTAL X (1.83) MULTIPLE SEQUENCE ALIGNMENT 

File: C:Documents and SettmgsK.V.BhatMy DocumentsKoIeKOLE SEQKSUSYD.ps Date: Thu Jun 04 10:20:20 2009 
Page 1 of 2 


halnl«n*_BB07__0lA 
hainian«JBBDl5 0 1A 
Biungo_BBOS OS 
ullvaftrli BB2541 
ailvaatria BBI.SJ 2* 
Bilv«*fcrI*_BB2£42 
mungo PDV1 
Bilvastri8_BBL50_2I 
«ongo_YMto*»2 
«ung6 TMV1 
■ungo IPB99 10 
lv*«trli_BlL7« 2000 
■ungo BBff3 10 
ailvaatria BBLSS 21 
jwngo_B8Di4_01B 
mungo BB414 
anmgo_BBOff8 01B 


halniaoa BB07 OU 
hainl*n»_BB0lS“O U 
mungo BBOJ Of 
#ilv*»tri« BB2?41 
silveatrla b5i» 83 2K 
silv«»fcrI*_BB2§42 
' mungo PD01 

•livMtxi* mkHo m 
wmgSjftMbini 
. mmgo flm 
*ungo tWJt 10 
Lvttatri* B8L74 2500 
*uigo BJ&33 10 
*ilv**trl* BBLSS"*'!* 
mmgo^BBW4 ffl» 
, mungo 00414 
mungo BBD5 s OU 

Illvwfcxii jmw® 2* 

mungo PBOl&K 
illvmatxi*' SBL40 2* 
<ilva«trii“lC2Sli(n 
halnlaiui BB2S4S 













1 1 2 1 25;' 

g$ nmm 

! mw 

I ini 

g g, hs 
* * $ * « a 

rf TO 7/TO 

,> | . • ■; ■ g 

? § ~ i : ' 3 ' 

'? <\ 7 r a - 7j \- , 

4 £ r i f c, ‘ 

» § ^ Jg ^ 

II mn 


' V x -fa,' t TO - TO pv TO TO 

vgjfS ffiM 

MSS' 


ippliP## 


s r ^ ? !" ■ ; fir* ■ ■ - * ■ 

M ff»» / JWT 

* ,‘fn^ *| 

to- *• 

’. 1 ? 7 ■ , 


itrry | gw wm* 
wr i? gj ^^fww 
fffff* ft gW TOWW 1 
Tfff ft ft»r* WfVfW 1 
yyfyr ft gw 
ffWff ft ftfW WPWW 
iw | gw * sm rt 
rr® % few «wfp 

«pfw ft TO JWWI: 

tw ft TO wrwr 
1W ft TO .WWW 
WV ft TO 
tfft® % to jwjtw 
ft TO *wm < 


„ , 3 - ^ 'x * i t * vf in? s? - „ r : ? v r^ f ; 

ft& ! ^ u ; f? *? to , , rfe V * * r : 

■zk t dft ? ft TO £ r t 7 u ** TO* : 
Jb | J|| ft ft TO ; J? P?fft \ w ft g rTOj 4 

■§* | 5^ | $? fefj ^ f| 4$? f > fjfft | 1 

s- SSTOTOto 1?te^M^TOTO ft rt ■■ sift ; : 

^■'TOjtoTOTO ^ Nro § ? : ;TO 

— ' ft ft®? ff ft TO o ‘/J ‘V? =f?TO ; V to f r; TOt f« 

to j Xp - ?J to ^ to ftt to ;. * „ 1 * !* V ^ r <V„ ft 1 

% P a/> »“ ^ ^ } / y~J % \f i u , i s~, t f e l - 

7k sgft f; 5.4 f i ' K 

*sT r? ts^ ft jf Yi ; V ^ ^ : B . ' Jr k 4 I#; 


S^S| v~,t ;«w ,i U« 

■JHt WfWt,***'! ,4 -v ^ ^ •■ f ^- 

CT 5 ftffi & f "^TO H 'j ’/* ’ ^ H' J " 

WM!| M|^«ivh^m n»| w V ' m 5% Wfc 

w*i^ ' -*■;•' * 

M0fj ^ t' 1 , > X 4 ml 

*” 5 ■ - t+s 

ftVr’TO Wj ^ r ,w ’ ,i! 

. . «;;.: 

i«4^ *»» • » h *v- 

mp Hxpap *,, ^ % l * % * «h 

"] l 1 ' ,u, TO 1 *' ‘ , ’' 

M/M * »> 

%-f? i b* 4 **Pjr*tf> »-V ( *fl 


-; ■: •■,■•':. '- v .v-. -';' v ';•: a ■' '■' 7 ^ ■ ,• ,/"}■■■ \ .•'■[ "■ ‘'■ 1 '* V ■ 

• v-.-: 

». 3 > rn - > ; ; " . : f- -r 4 

»* ..■ .,. i-:s> .. '.;f- to to 

: v.'rr : ;7'V' : , ■ 4 :aV."-' r.v-; ;, : ' 

to.toto-.:- : ,.: ' .... v..:. ;., U v • ■ ■ 

.~4 to- YY; 3 r: -■■■o- : % 7 : 

>> * i. h ;* , N» #*■ 

- to ''\;''7' ., '• ' ; ; - v. ;r .;. ■'■.;■ “. -• Y, V"' 

-;.-- ■ , ' .'■ .'■ '" ' ■. - ■ •.. - : :: -. ■ ' - • • • - ' : ' '^,to ' 7 ' ' Y '■' " 

' ?? ' toTO' fT %. \ .c^- • . ^v"- 

•>“•■ to-fi; :• '• '• ; ^ ; r ■-' - toto :■’. ;; >7 .7'- 7: 


:<?AGMCQM<3AACA< 


TOSwS 

mungo_i 

6 ilvoctri^BBL5| 
«ungo_Varei 
arungoto 
*ungo XPDS! 

w AW, 

a uiigo_BBO. 

,ilY*tri* BBL5! 

jHingo BBB14. 

mungo_u< 
mmga BBD08 

BBL4( 

;SIStri.:i c«] 

baiiiix«na_BBi 


•iniana BB07 01A ACCOCXCCflftOCC 
S— ACCGC^GGMOCC 
JSSioW 09 ACGOC|COfcMCC 
n^SSii BB2641 AGCOC CCUUMCC 
2”S bBL 63 2K ACCflC CflWSCC 
nJ«trrf BB264.2 ACCOC CGAAOCC 
‘ wo PDTJ1 ACC«C CGttflCC 
BBLSO 2K ACCOC CCttAGCC 
Ss.Vambina ACCGC CGMGCC 
nungo TMVl ACCOC COAAOCC 
m ungo IPUfd 10 ACC(rfGOW<JCC 
,t5I BBL76 2000 ACCGC COAAOCC 
»uHgo BBS* 10 ACCGC COAAOCC 
,„?ri ; BBL55-2K ACCOC COAAOCC 
muigo BBD14 01B ACCOC COAAOCC 
muSo 00414 ACCOC COAAOCC 
mungo BBD08 01B ACCO^COAAOCC 
BBL80 2K ACCOC COAAOCC 
mungo PDU1016 ACCOC COAAOCC 
r«?rii SBL40 2K ACCO^COAAOCC 
706 trickle 25 14 07 ACCOCACOAAOCC 
hainiaHa_BB2649 ACCOCACOAAOCC 
rul«c 460... 


CLUSTAL X (1.83) MULTIPLE SEQUENCE ALIGNMENT 

File: C:Documents and SettingsKV.BhatMy DocumentsKoIeKOLE SEQKSUMTD.pDate: Thu Jun 04 10:18:57 2009 
Page 1 of 2 


nungo BBB8 01B - 

mungo_IP09 9_1 0 - 

lvnafcr i a_BBL4 0_2K - 
mungo_TMVl - 
mungo_Vanib«n2 - 
®ungo_BB03_09 - 
nungo PDU1 - 
lv*»tri*_IC2E>140? - 
nungo_OQ4l4 - 
lv»»fcri«_BBL63_2K - 
8ilv'*itrl8_BB2642 * 
Ivoabr i »_BBL8 0_2K - 
lwstrie_BBI.50 2JC - 
oilr*etri» BB264X <J? 
lv«8tri8_BBIi7 6_2K - 
rulor 1 


nungo_BBB8 01B 
iaungo_IP095_l 0 
lv**tri87BBL40_2K 
munga_THVX 
»ungo_Va5b«n2 
nungo7 BB 08_09 
mungo PDU1 
lv*8bri8_lC25l407 
nungo_D0414 
lv«»tris BBL63 2K 
8ilVft8trIfl_BB2?42 
Xv*8brl«_BBLe0 2JC 
Xv«8tri8_BBIi50_2K 
silr«0tris_BB2641 
Xv*stri8_BBL76 2K 
rulnr 


nungo BED 8 01B 
nungo7lP09S5_10 
Ivesbrl s3bbL 4 032K 
mungo_TM91 
mungo_Vaiban2 
mungo_BB03 09 
nungo PDOl 
lv»»tri *_IC2§1407 
nungo UG414 
lv#8tri8 BBL63 2K 
ailvnsbrl* BB2642 
lvnstri 8 _BBXi 80 2K 
lv»*tri» BBL50~2K 
ailvnitrls BB2641 
1 v*b tr i «_BBL7 6 S9C 
' rulor 





E Twiwm 

ff 8§ «; 5*fl rtm few 5 

SSfSUSRSSS |5 
§& 

illtl ;3 1 ■ ff 7* Iff 

gg^rJ^ n gs ? ggj te 

WOTW& *< H « ?p 

g^fs4-wfe "i % * gyr 


■ v-.'- SW JWffiS 
• ipp ;gwr ,tn», 

■’ fffff*f?V 

7 CTfTF! Wf ^W#p. 

■ -y?- ■-’ ;/f ’ “jWfft*? 


CLUSTALX (1.83) MULTIPLE SEQUENCE ALIGNMENT 


File: C:Documents and SettingsKV.BhatMy DocumentsKoleKOLE SEQKSHMTD.pDate: Thu Jun 04 10:18:57 2009 

Page 2 of 2 


wngoJBBD8 01B 

nungo IWM_10 

;S!OBL40„2K 
nungo_ TMV7. 
nungo Vamb*n2 
mungoJ3B0 3_0 9 
nungo PDOl 
.CrL-KSfllO? 

nungo_u«4X4 

; tric 8BL63„2K 
r *0fcrIa„BB2«42 
( tria_BBL80_2K 
(tri»~BBLS0 2X 

itrin BBIi7 6_2K 
~ rulor 


rnungo BBD8 01B 

«ungo~iPD99_io cp:c|COTa«KCCc®k 
lv«»tri*_BBL40_2K ^ 

mungo_TMVl ^yCy^|Q^CCC<3- - 

mungo _V«mb*n 2 

mungo_BB03_09 
mungoJ0G4l4 

liiltei?BBLOO_ 2 E - - - - - - --- ------ 

1 vo ctrie BBL50 2K 

Svwtri._BB 2 ? 4 i mmmmmm 

lvostrio BBL7 6_2K ----------- --- 




•CCC«M88- 



discussion 



Advances in cellular and molecular biolosv have , , ... 

gy nave provided new tools like marker aided 

selection, which increases the efficiency of breeding programmes especially disease 
resistance breeding. In conventional methods, selection of genotype relied on field or 
glass house screening. It may also require the use of controlled environment, when the 
disease development is temperature sensitive as that of yellow mosic virus. Hence, the 


screening process is complex and time consuming. But in molecular marker assisted 
selection, once a gene of interest is tagged with molecular marker, selection for that gene 
can be made based on the marker. In present study, an attempt was made to identify the 
RAPD marker associated with MYMV resistance, using bulk segregant analysis. One of 
the most time consume requirement in marker development is the need to screen entire 
mapping population, suggest that testing the entire population is required only when 
polymorphism between the bulks are detected. This results in a considerable saving of 
time particularly when used PCR based technique such as RAPD (Williams et al. 1990). 

A total of sixty six random decamers were surveyed for identification of 
polymorphic markers between the parents and twenty four produced specific band for 
resistance parent which were absent in susceptible parent. However the parents were 
chosen at Allahabad according to there specific characters such as resistance for MYMV 
and it were clearly distinguished there. But at Delhi all the F 2 individuals were found to 
be susceptible along with the parents except some individuals (Plant no 47, 65, 80 and 90) 
according to interaction phenotypic scale. Due to some unavoidable conditions and the 
total number of F 2 individuals was found to be very less. The climatic conditions in which 
the F 2 individuals were grown was not appropriate. Either due to pathogens found in Delhi 
or the environmental conditions the number of F2 was reduced to 46. These small number 
of F2 population were not sufficient to generate a molecular map. Hence our study was 
further extended to find the closest wild species related to our crop which are resistance to 
some disease so that diversity analysis and gene flow can be studied. Further we tried to 
search for the single nucleotides which made differences. 
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W,ld relatives of crops species are source of important genes for agriculture ^ 
genetic divereity measured a. the biochemical and DNA level is generally greater in wild 
species titan there related cultigen (Hallan, 1984; Xu e, al„ 2000). Intensive modem 
breeding has contributed to a narrowing of crep genepools as a few improved cultivere 
dominate large areas ( Ladizinsky, 1985). Due to recent development of gene transfer 
technology, genes from cross incompatible wild species can be use for breeding (von 
Bothmer & Seberg, 1995). Therefore, evaluation of a broad array of wild species is an 
important approach to seeking genes which rare or absent in cultigens. 


Wild Relatives as Genetic Resources 

Fuzi and Miyazaki (1987) reported an accession (TC1966) of V. radiata van sublobata 
that showed perfect resistance against azukibean weevil ( Callasabruchus chinensis). The 
resistance was found to control by a single dominant (Kitamura et al, 1988). Fuzi et al 
(1989) further found that TC1966 is completely resistant against C. maculates, C. 
phaseoli and Zabrotes subfasciatus. Tomooka et al (1992) developed a bruchid resistant 
mungbean line in Thailand by using TC1966 as a gene source. In addition to bruchid 
resistance, high resistance to yellow mosaic virus (Singh and Ahuja, 1977), high 
methionine content in seeds (Babu et al, 1988), higher photosynthetic activity and 
tolerance to drought ((Ignacimuthu and Babu, 1987), higher tolerance to saline and 
alkaline soils (Lawn et al, 1988) have been reported for V. radiata var. sublobata. In 
comparison to V. radiata var. sublobata there are very few studies of evaluation of V. 
mungo var. silvestris as genetic resource. However, this variety is cross compatible with 
V. mungo (Miyazaki et al, 1984). Therefore, further collection and evaluation V. mungo 
var. silvestris is necessary for genetic upgradation of V. mungo. 

Information on the levels and distribution of genetic diversity of any plant species 
may contribute to the knowledge of their evolutionary history and potential, and is critical 
to their conservation and management (Hamrick and Godt, 1996). Several markers may 
be used to identify and assess the genetic diversity and phylogenetic relationship in plant 
genetic resources. The traditional method based on morphological traits requires 
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extensive observation of mature plants and can 

n not serve as unambiguous markers 

because of environmental influences (Wrigley et al , 1987) 

With the advent of molecular markers ^ non, 

? a new generation of markers has been 

introduced over the last two decades, which has revolutionized the entire scenario of 
biological sciences. DNA-based molecular markers have acted as versatile tools and have 
found their own position in various fields like taxonomy, physiology, embryology, 
genetic engmeenng, etc. They are no longer looked upon as simple DNA fingerprinting 
markers in variability studies or as mere forensic tools. Ever since their development, they 
are constantly being modified to enhance their utility and to bring about automation in the 
process of genome analysis. The discovery of PCR (polymerase chain reaction) was a 
landmark m this effort and proved to be a unique process that brought about a new class 
of DNA profiling markers. This facilitated the development of marker-based gene tags, 
map-based cloning of agronomically important genes, variability studies, phylogenetic 
analysis, synteny mapping, marker-assisted selection of desirable genotypes, etc. Thus 
giving new dimensions to concerted efforts of breeding and marker-aided selection that 
can reduce the time span of developing new and better varieties and will make the dream 
of super varieties come true. These DNA markers offer several advantages over 
traditional phenotypic markers, as they provide data that can be analysed objectively. In 
this article, DNA markers developed during the last two decades of molecular biology 
research and utilized for various applications in the area of plant genome analysis are 
reviewed. 

A total of 375 individuals (15 seeds of 26 accessions of 3 species) of vigna 
showed the polymorphism within accessions and within groups using STMS analysis. The 


ten polymorphic primer pair used to study in which two loci were found to be highly 
polymorphic with 5-4 alleles. V. hainiana and V. silvestris both showed maximum 
number of polymorphism. The product size ranged from 150bp to 320 bp. This showed 
the efficiency of the marker used and the fingerprinting methods. However, compared 
with the use of morphological markers, the results obtained using microsatellite markers 
are more stable because the repeat sequences are hardly affected by factors such as the 
environment, sampling time, which part is sampled and especially the choice of external 
factors during evolution. 
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Microsatellite analysis needs little DNA hut ,. i 

NA but ’ detects much polymorphism, gives 

reliable results and is simple to operate. Therefore, STMS is a convenient and reliable tool 
,o quickly identify the genetic diversity of large numbers of gennplasm (Hong- Liang et 
al„ 2004). Specific markers like STMS (sequence-tagged microsatellite markers) or STS 
markers have proved to be extremely valuable in the analysis of gene pool variation of 
crops during the process of cultivar development, and classification of germplasm. These 
markers are extremely sensitive and can detect allelic variability during cultivar 
development. STS markers specific to chloroplast or mitochondrial DNA have been 
useful in providing seed and pollen specific markers which can be utilized for the 
detection of length variation at multiple physically linked sites and may be used to 
provide haplotype data and thus genotypically unique individual plants. 

Population structure is the consequence of multiple events, starting from the 
modalities of domestication, the structure being further shaped by modem breeding. The 
lack of clear population differentiation, high heterozygosity and gene flow, low fixation 
indices all point towards a high rate of germplasm migration across geographic regions in 
the country. There is a considerable exchange via buying, sharing of cultivars across 
mungo growing belts as well as introduction and spread throughout the country. Analysis 
of molecular variance, low FST value, cluster analysis and high values for estimates of 
Nm indicated presence of substantial gene flow between the cultivars from different 
regions, possibly due to cultivar exchange followed by gene introgression by crossing and 
selection. 

The POPGENE analysis was done for the populations as well as for the groups. 
Groups were classified into three. Group I was formed of V. mungo. Group II was of V. 
silvestris and group III was of V. hainiana. Total number of alleles, effective number of 
alleles and Shannon’s Index were calculated for 25 accessions as well as for three groups. 
The effective number of alleles varied from 3.0 to 1.7. Shannon’s index was maximum in 
silvestris BB2641. In three groups Shannon’s index was highest in V. m. silvestris which 
is a wild relative of Vigna mungo. 

The observed heterozygosity and expected heterozygosity were also calculated for 
population and groups. The mean values of observed and expected heterozygosity for 25 
populations were 0.084 and 0.335 respectively. Since observed was less than expected the 
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values suggested that there was no cross pollination also with ; , 

* on also Wlth m populations. For groups 

the observed was less than expected suggesting self pollination nature of tndbean. 

Nei’s original measure of genetic distance was calculated for population as well as 
for groups. In group I the closest relation was found in tnungo BBD08-1B, IPU99/10 
TMV 1 and PDU1 . For group II the closest relationship was found in between V. stives.* 
BBL63-2Kand BBL80-2K. IN group III the accession BB07-2k and BBD15-2K was the 
closest. Population dendrogmn was made by Nei’s genetic distance. All 25 accessions 
were clustered. Dendrogram were clearly clustered no inter groping was found in Ms 
dendrogram. Dendrogram were clearly clustered no inter groping was found in this 
dendrogram. Two main clustets were formed with three species. Ouster I consist of V. 
tnungo and V. silvestris accessions and cluster II was V. hainiana accessions. Cluster 1 
again divide into two, on which one contain V. mungo with some accessions of V. 
silvestris. Another one consists of four V. r. silvestris accessions which were closely 
relative. 


Fixation statistics were produced for individual SSRs and groups of germplasm. 
By using 1000 permutations significance of the estimates were obtained. The percentage 
of variation within populations (16.03%) followed by 39.13% among populations within 
group and 44.84% among groups. Average fixation indices were FSC: 0.710, FST: 0.840 
and FCT. -0.84 indicative of high level of outcrossing and gene flow and no geographic 
isolation at any hierarchical level. High levels of gene flow via germplasm exchange may 
have caused orchard stands in the northern plains and related regions to be genetically 
connected with little genetic differentiation. 

The population differentiation value (Fst) was used to determine the amount of 
geneflow (Nm). The values were calculated separately for populations and for groups. Fst 
value is inversely proportional to Nm. The Nm value for the populations was found to be 
negligible, suggesting no gene flow within populations. For the groups, Nm value for V. 
mungo was 0.4647 which was more than both V. r. silvestris and V. hainiana. Since V. 
mungo is a cultivated it showed the highest gene flow within. The Ewen-Watterson test 
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for neutrality was done using 1000 simulated samples for the in w „ 

pics I0r tne 10 loci all were neutral to 

selection pressure. 

The numbers of segregating s«es (S), observed nueieotide dive^ty per site 

between any two sequences assuming tha, the sample is tandom (p), nmnbe, of 

haptotypes (H). haplotype diversity (Hd), and average number of pairwise nucleotide 

diffetences within population (K) were calculated using DnaSP. Tajima's D expend ,o 

positive natural selection whereas negative value* . 

f negative values correspond to negative or purifying 

selection. The genetic differentiation among the parasite populations was calculated in 

terns of fixation index (Fs.) that estimates divetsity within a subpopulation with respect 
to total genetic 

diversity. In addition, average number of pairwise nucleotide differences (Kxy), 

nucleotide substitution per site (Dxy), and net nucleotide substitution per site (Da) 

between populations were also calculated. The above parameter* were also 
estimated on DnaSP. Phylogenetic analysis was perfotmed by neighbour-joining (NJ) 
method with Kimura 2-parameter distance matrix inMEGA version 3.0 

For the sequencing of the samples sequence tag microsatellites (STS) marker was 
used. Mostly haplotype diversity, Tajima’sD observed pairwise nucleotide diversity was 
calculated. Positive value for locus OLM1 of tajima’s D indicates the positive neutral 
selection. Measure of DNA sequence polymorphism and neutrality test table(table 16) 
indicates that OLM1 marker was more informative than other markers. Fst values for 
Loci SUSY was highest 0.628 and Nm was calculated to be maximum 7.11 in the locus 
OLM1 indicated that the geneflow was very high mainly suggesting some cross 
pollination in the crop. 

Inter population nucleotide difference (Kxy) varied from 0.0235 (between 
silvestris and mungo) in locus UDPGDB3 to 26.639 (between mungo and hainiana ) in 
OLM1. Similarly, the average number of nucleotide (Dxy) and net nucleotide (Da) 
substitutions per site between populations ranged from 0.00088 (SHMT1) to 0.09146 
(OLM 1 ) respectively. Inter-population comparison of the total (TM) versus shared (SM) 
number of mutation showed only on locus OLM1 (29 of 53 mutations were shared) 
(Table 17). Gene flow was varied from 0.31 (SUSY8) to 8.83 (OLM1) on basis of 
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haplotype data —on (Table ,8). Based „„ sequence data InfoCon i, was found 
0.15 (SUSY8) to 7.1 1 (OLM1). Total number of segregating sites (S) for Population wise 
Measures of DNA sequence polymorphisms and tests of neutrality was varied fem 16 

(RNAR 8) to 5 (UDPGDB3). Number of haplotypes also high 5 in RNAR 8 where for the 
r^t two locus were 3 (UDPGDB3) and 4 (SUSY8) respectively. Haplotypes divemity and 
nucleotide divemity value were 0.53216 and 0.01299 which is very high value than otirer 
two loci. Individually vigna hainiana showed maximum value in all respect for all the 
three loci. Above information showed the relation between cultivate and wild species. In 
the all aspect mungo and hainiana were difTer from each other. So for the future studies 
wild species can be used for breeding and crop improve programme. 


The gene trees were prepared for the locus SHMT, SUSY8, BV16 UDPGDB3 and 
OLM1. Three types of trees were discussed, Neighbour Joining, Minimum Evolution and 
Maximum Parsimony. For the locus SHMT, SUSY8 BV16 and UDPGDB3 neighbour 
joining tree and minimum evolution tree were almost similar in the pattern of clustering. 
Except for the arrangement of some species both the trees showed no variation in the 
cluster formation. Maximum parsimony tree showed the variation in the single nucleotide. 
Instead of ‘T’ either ‘A’ or ‘G’ was inserted. Similarly for the locus SUSY and SBV16, 
the cluster formation for neighbour joining and minimum evolution was similar. The 
maximum informative locus was OLM1. 

The alignment of the species for the three loci showed the single nucleotide 
insertions and deletions. C luster formation was done on the basis of this single nucleotide 
polymorphism. Very less number of mutations was observed hence variation was very 
less. Comparing between five loci OLM1 showed more information than SHMT, 
UDPGDB3, SUSY and SBV16 hence these locus gave more information. 
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SUMMARY 


Pulses are unique crops that are a rich source of wo^m 

ce ot vegetable protein and have in-built 

mechanism to fix atmospheric nitrogen. India is a maior 

6 a ma J° r P^se growing country in the 

world, with about 37 /o of area and 28% of the product™. Besides tte immense 
importance with regards to nutritional and food security, pulses have been uadMoually 
grown under margmai conditions of poor soil moisture, fertility and different other 
abiotic stresses, while the more of the fertile lauds, as well as the improved agmnomic 
practices have always been reserved for the cereal crops. These constraints have resulted 
i„ lower productivity and availability of pulses, besides significant efforts for improving 
the cultivais for fertile lands. One of the approaches to increase grain legume production 
is to expand the cultivation of different species in the advemely affected soils through 
developing cultivars which are more tolerant to different abiotic and biotic stresses. 


Mapping populations are usually obtained from controlled crosses. Decisions on 
selection of parents and mating design for development of mapping population and the 
type of markers used depend upon the objectives of experiments, availability of markers 
and the molecular map. The parents of mapping populations must have sufficient 
variation for the traits of interest at both the DNA sequence and the phenotype level. The 
variation at DNA level is essential to trace the recombination events. The more DNA 
sequence variation exists, the easier it is to find polymorphic informative makers. When 
the objective is to search for genes controlling a particular trait, genetic variation of trait 
between parents is important. If the parents are greatly different at phenotypic level for a 
trait, there is a reasonable chance that genetic variation exists between the parents, 
although uncontrolled environmental effects could create large phenotypic variation 
without any genetic basis for the effects. 


Blackgram or urdbean is one of the important grain legumes of India and 
Munhgbcan yellow mosaic virus (MYMV) disease is one of the most devastating 
diseases. Munhgbcan yellow mosaic virus (MYMV), a whitefly-transmitted 
bigcminvirus, causes disease in blackgram. This virus is assigned to the genus 
Begomovirus within the family Geminiviridae. Plant resistance is the most common and 
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efficient method of MYMV control. In blackeram w 

, , . ’ symptom types yellow mosic and 

necrotic mottle, can be distinguished (Nair and Nene iq™ A 

. f ’ ^4). A strategy to improve the 

efficiency of mapping, named selective manning 

PP g, was proposed by Visionet al. (2000). It 
consists of a two-step process in which, first a mannw , - 

’ a ma PPmg population of usual size (Ny 4 60- 
250) is used to construct a saturated framework 

P Wl ^b markers placed on it with high 
precision, and second, new markers are added to tw* 

, , 6Cltothlsma P with lower precision using a 

selected subset of highly informative plants. 

To attain the objectives of mapping MYMV resistance with DNA markers, i, was 
proposed to deveiop F 2 popuiation from a cross between IPU 982, an improved high 

yielding recommended variety but highly susceptibie MYMV infection, and Cudaclc 

Local, a local landraces (LLRs) from the Cuttack c r, ■ 

) ic v^unack district of Onssa known to be a low 

yielder but resistant to MYMV disease. The erne* a 

J ne cross was done at Allahabad and the FI 

progenies were also planted there. The F2 general,-^ • 

F generation was grown m NBPGR, New Delhi 

along with the parents. But at Delhi all the F 2 individuals were found to be susceptible 

along with the parents except some individuals (Plan, n0 47 , 65> 80 md %) accorifag , 0 

interaction phenotypic scale. Due to some unavoidable conditions and the total nmnber of 
Ft individuals was found to be vtny less. The climatic conditions in which the Ft 
individuals were grown was not appropriate. Either due to pathogens found in Delhi or 
the environmental conditions the number of F2 was reduced to 46. These small number 
of F2 population were not sufficient to generate a molecular map. Hence our study was 
further extended to find the closest wild species related to our crop which are resistance 
to some disease so that diversity analysis and gene flow can be studied. Further we tried 
to search for the single nucleotides which made differences. Some more accessions were 
used and the main objectives of using these markers were 


v) Elucidation of genetics of resistance to YMV diseases in urdbean 

vi) Molecular mapping of genes controlling resistance to YMV. 

vii) Studying genetic relationships among important Vigna cultivars and wild 
species using molecular markers. 
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iv) Study of genetic variability in different oe n( => u 4 

t gene pools based on molecular 

gene sequence 

DNA for this analysis extracted following the CTAB exaction protocol of Saghai- 
Maroof a al (1984).The STMS technique was optimised for DNA quantity (40 ng) 
Mg a 2 and annealing temperatures of the primer pairs. Initially 66 RAPD prime* were 
chosen but later on the more efficient makers were used. A total of 10 STMS and 6 STS 
primer pairs generating optimum and consistent PCR profiles were selected for the 
analyses. The DNA profiling of the 25 accessions with selected primer pairs yielded good 
polymorphic amplification products. The PCR preducts that were not consistent with the 
estimated allele sizes for a locus were not scored since these were considered to be non- 
specific amplifications. 


. The data was also used to analyze the population genetic parameters such as 
total number of alleles, effective number of alleles, Shannon Information index etc. The 
maximum number of alleles found was 5 with V. r. silvestris. Ewens -Watterson test for 
neutrality was done to analyze the occurrence of selection and gene flow. Analysis of 
molecular variation (AMOVA) was done to find out the genetic differentiation between 
and within populations and among groups. Dendogram was made for the 25 accessions to 
show the arrangement of the accessions cluster wise. 


The sequencing was done for the samples by which the numbers of segregating 
sites (S), observed nucleotide diversity per site between any two sequences assuming that 
the sample is random (p), number of haplotypes (H), haplotype diversity (Hd), and 
average number of pairwise nucleotide differences within population (K) were calculated 
using DnaSP. Tajima’s D correspond to positive natural selection whereas negative 
values correspond to negative or purifying selection. The genetic differentiation among 
the parasite papulations was calculated in terms of fixation index (Fst) that estimates 
diversity within a subpopulation with respect to total genetic diversity. In addition, 
average number of pairwise nucleotide differences (Kxy), nucleotide substitution per site 
(Dxy), and net nucleotide substitution per site (Da) between populations were also 
calculated. The above parameters were also estimated on DnaSP. Phylogenetic analysis 
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was performed by neighbour-joining (NJ) method with Kimura 
matrix inMEGA version 3.0. 


2-parameter distance 


The major conclusions drawn were tha, due recent development of gene 
transfer technology, genes from cross incompatible wild species can be use for breeding 
Therefore, evaluation of abroad anuy of wild species is an important approach to seeking 
genes which rare or absent in cultigens. 
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ANNEXURE 



annexure-i 

1. Solutions, chemicals and reagents used for DNA extraction 

1. Liquid Nitrogen 

2. Cetyl T rimethyl Ammonium Bromide (CTAB) Buffer for DNA eittraetiou: 

(i) CTAB (10%) 

1 Og CTAB was dissolved in sterile distilled H 2 0 and volume was made upto 100 

ml with distilled water. 

(ii) Sodium chloride (NaCl, 4M) 

292. 2g of Na( l was dissolved in distilled H 2 O and volume was made upto 100 

ml. The solution was autoclaved prior to use. 

(iii) Tris: Cl buffer (pH 8.0, 1M) 

1 2. 1 1 g of Tris salt was dissolved in distilled H 2 0 and volume was made upto 100 
ml and pH was adjusted to 8.0 using 1 N HC1. The solution was autoclaved prior to use. 

(iv) Ethylene Diamine Tetra Acetic acid (EDTA, 0.5 M) 

18.62 g EDTA was dissolved in sterile distilled H 2 0. The pH of the solution was 
adjusted to 8.0 using IN NaOH. The volume was made upto 100 ml using sterile distilled 
H ; C) and the solution was autoclaved. 

(%) 2-Mercaptoethanol (2%) 

2 % solution provided by manufacturer was used directly. 


Buffer composition 




Component 

Stock sol. 

Working buffer 

Vol. of stock taken to 

prepare 200 ml buffer 

UTAH 

10% 

1.5% 

40 ml 

2. NaCl 

4M 

1.4 M 

70 ml 

3. I ns 

INI 

lOOmM 

20 ml 

4. EDTA 

0.5 M 

20 mM 

8 ml 

5 Mercaptocthanol 

2% 

2% 

4 ml 

6 Distilled H O 

- 

- 

58 ml 


Total 


200 ml 


3. Isopropanol 

4. Sodium Acetate solution (3M, pH 5.6) 

30.75 g sodium acetate was dissolved in sterile distilled H 2 0, pH was adjusted to 

5.6 with glacial acetic acid and volume made upto 50 ml. The solution was 
autoclaved and stored till use. 

5. Chloroform: Isoamyl alcohol (24:1) mixture 

96 ml of chlorofom, was mixed with 4 ml of isoamyl alcohol. It was stored in 
amber coloured bottle. 

6. 70% Ethanol 

70 ml of absolute ethanol was mixed well with 30ml of sterile water and stored in 
a stoppered bottle till use. 

II. DNA purification 

1. Phenol: Chloroform: isoamyl alcohol (25:24:1) mixture 

1 00 ml of Tris saturated phenol was added to a mixture of 96 ml chloroform and 4 
ml isoamyl alcohol. The mixture was mixed well prior to use and stored in amber 

coloured bottle. 

2. RNaseA (20 mg/ml) solution 

RNaseA : 20 mg 

Tris-Cl (pH 7.5) : 10mM 

NaC'l : 15mM 

Sterile water was added to make the volume to 1 ml. The solution was heated at 
1 (X)“C for 15 minutes to inactivate any DNase present and then stored in aliquotes at 
2()T. 

3 . Pronase ( 1 0 mg/ml ) solution prepared in sterile distilled water. 

III. Solvent for DNA 

I ris: EDTA (TF.) butler (10 mM Tris: ImM EDTA, pH 8.0) 

10 ml of Tris ( 1 M) buffer, pH 8.0 and 0.2 ml of 0.5 M EDTA, pH 8.0 was mixed 
with sterile distilled H 2 0 and volume made upto 100 ml. The solution was autoclaved 

prior to use. 
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IV DNA Quantification 

1. Hoechst Dye (H-33258) 10 X solution: 

10 mg of Hoechst 33258 dye was dissolved in sterile distilled H 2 0 and volume made 
upto 1 00 ml and stored in an amber coloured bottle at 4°C. 

2. 10 X TNE buffer stock solution (100 mM Tris: lOmM EDTA: 2M NaCl, pH 7.4) 
12.1 1 g Tris, 3.72 g EDTA and 116.89 g of NaCl were dissolved in sterile distilled 

water and volume was made upto 1 00 ml using distilled water. The pH was adjusted to 
7.4 with cone. HC1 solution, filtered before use and stored at 4°C. 

3. Assay solution 


Component 

Low range DNA assay 

High range DNA 


10-500 ng/ml final DNA cone. 

assay 100-5000 ng/ml 

Hoechst 33258 stock soln. 

10 pi 

100 pi 

10 X TNE 

10 ml 

10 ml 

Distilled H;0 

90 ml 

90 ml 


prepared fresh each time 
4. DNA standard 


i. Low range assay: 

1 pl/ml of calf thymus DNA standard was used at 1:10 dilution (100 pg/ml). 2ml 
of calf thymus DNA (100 gg/ml) was mixed with 2ml assay solution for low range which 

gave 100 ng/ml standard solution. 

ii. High range assay 

2pl calf thymus DNA standard (1 pg/ml) was mixed in 2 ml assay solution for 

high range assay gave 1000 ng/ml standard solution. 

V. GEL ELECTROPHORESIS 
1 . Agarose gel ( 1 .8%) 

4.5 g agarose was added to 250 ml with 1 X TAE buffer, the contents were mixed 
thoroughly and boded for 2-5 minutes to dissolve the contents. The mixture was cooled 
down to 4UT The molten gel was cast in a gel tray with a comb containing 33 teeth to 

produce wells. 
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2. Ethidium bromide (10 mg ml' 1 ) 

10 mg of ethidium bromide was dissolved in sterile water and volume made up to 1 
ml- The solution was stored in an amber coloured bottle, at 4°C 

3. Loading dye (1 OX) solution 


1 . Bromophenol Blue 

2. Xylene cyanol FF 

3. Glycerol 

4. TAE 


0.25% 

0.25% 

50% 

IX 


4. Tris: Acetate: EDTA (TAE) buffer - 50 X (stock) solution. PH 8.0 

2M tris-acetate, pH 8.0 
0.05M EDTA, pH 8.0 
VI. PCR COCKTAIL 

1. Taq DNA Polymerase 

A stock solution of 3 units pi' 1 was provided by the manufacturer (Bangalore Genei) 

was stored at ~20°C. 

2. 10 X Assay buffer 

10 X PCR assay buffer for Taq DNA polymerase containing 15mM pi' 1 
magnesium chloride provided by the manufacturer (Bangalore Genei) was used. Storage 

was at 20T. 

3. Deoxy ribonucleotide T riose phosphate 

dATP (10 mM). dGTP (10 mM), dCTP (10 mM) and dTTP (10 mM) were mixed in 

equal volumes and stored at -2(fC till use. 

4. Magnesium chloride 

A solution of 1 5 mM pi 1 provided by the manufacturer, stored at -20°C, was used. 

5. Primer 

The primer was provided by the manufacturer in a lyophilized form. Based on the 
molecular weight of a given primer, a solution of 6 pM was prepared by adding the 
required amount of sterile water. Storage was at -20 C. 


IV 


Chemicals and Materials used 


Type 

Material 

Source 

Molecular 

weight 

markers 

lOObp DNA ladder 

1Kb DNA ladder 

Fermentas 

Fermentas 

Modifying 

enzymes 

Taq DNA Polymerase 

RNase 

Fermentas, Applied Biosystematics 
NEB 

Qiagen 

Kits used 

Dneasy DNA Extraction kit 

PCR product purification kit 

Agarose gel DNA purification kit 

Qiagen 

Bangalore Genei 

Himedia 

Dyes 

Ethidium Bromide 

Xylene cyanol 

Methylene blue 

Amersco 

Sigma 

Sigma 

Fine 

chemicals 

APS, EDTA, CTAB, Acrylamide, 

Bis- Acrylamide, TEMED, Agarose 

Sigma, Amersco 

Locally 

available 

chemicals 

Isopropanol, Iso-amyl alcohol, 

CaCl 2 , NaCl, NaOH, Glucose, 
Methanol, MgCI 2 , Chloroform, 
Glycerol, Acetic acid, MgSC> 4 , 
Formaldehyde, {3-Mercaptoethanol, 
Tris, Ethanol 

Qualigens, Sigma, SRL 
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annexure-ii 


Table 5. Allele frequency of 25 vigna populations. 
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Table Allele frequency of three groups of rigna. 
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